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ABSTRACT 

The objective of this work was to determine the expected 
deflections prior to and at the rupture load for a panel of hull plating 
of an ice-strengthened ship. An idealized mathematical model of a pin- 
nacle of ice being forced against the plating was used as the assumed 
loading mechanism. The panel of hull plating was modeled as a clamped 
rectangular plate. The idealized ice-loads were defined as a four sided 
pyramid of ice which when forced against the center of a rectangular 
plate crushed at a given pressure. Thus the size of the load bearing 
surface continually increased as the applied force was increased. 

A method of analysis was established which provides a step-by- 
step solution of the large deflections of a plate and a means of deter- 
mining the minimum expected failure load. The analysis is begun at the 
point where the load on the plate from the pyramid of ice is large enough 
to cause plastic collapse by a three hinge mechanism. A solution method 
for the deflection of a uniformly loaded membrane forms the basis for the 
calculations of the plate deflections. Then the energy absorbed by the 
plate at the point of expected failure is calculated from an equation 
that depends on the total strain at the point of failure. 

The possibility of failure of the plate at any given load and 
deflection may be determined from the failure criterion. This criterion 
was established by forming a relation between a parameter representing 
the geometric properties and a parameter representing the material pro- 
perties of the plate. The relation was formed from published data on 
plates that had been tested to failure. The calculated failure parameters 
for these plates were used to form a linear equation which correlates 
all the experimental failure loads. The criterion was then checked by 
comparing the results given in other experiments on the large deflections 
of plates. 



The method of analysis was then applied to plates similar to 
the hull plating of ice-strengthened ships. The results showed that the 
method of analysis can be used with a reasonable degree of accuracy to 
calculate the large deflections of a plate subject to a lateral load and 
to predict the load and deflection at which failure may be expected. 

Thesis Supervisor: J. Harvey Evans 

Title: Professor of Naval Architecture 



2 



ACKNOWLEDGMENTS 



I would like to express my appreciation for the guidance and 
encouragement given me by my supervisor, Professor J. Harvey Evans. I 
would also like to thank Associate Professor Norman Jones for his assis- 
tance and insight on the behavior of plates. 

The foremost acknowledgment is for the patience and understand- 
ing of my wife and family. 



3 



TABLE OF CONTENTS 



Page 

ABSTRACT 2 

INTRODUCTION 11 

Discussion of Problem 11 

Idealization of Ice-Loads 12 

Type of Analysis 12 

Load/Deflection Relations 13 

Failure Analysis . 14 

Outline of Procedure 15 

Ordering of the Bibliography 15 

PROCEDURE FOR CALCULATING ICE/PLATE INTERACTION 17 

Initial Assumptions for Ice-Loads 17 

Initial Assumptions for Rectangular Plates 17 

Model of Ice-Loads 20 

Plastic Collapse Load 21 

Plate Deflection Calculations 23 

Load /Deflection Relations 27 

Energy Absorbed Relations 29 

Failure Criterion 33 

Summary of Method of Analysis 37 

RESULTS OF APPLICATIONS OF THE METHOD OF ANALYSIS 40 

Plastic Collapse Load and Deflection . 40 

The Membrane Solution for Plate Deflections 41 

Energy Absorbed Relation 44 

Application of Failure Criterion 46 

DISCUSSION OF RESULTS 54 

Load/Deflection Relations 54 

Absorbed Energy Relation and Failure Criterion 59 

CONCLUSIONS 61 

RECOMMENDATIONS 61 

Appendix A 

PLASTIC COLLAPSE LOAD 62 

Appendix B 

DERIVATION OF MEMBRANE FORMULA FOR PLATE DEFLECTION... 74 



4 



TABLE OF CONTENTS (CONTINUED) 



Page 

Appendix C 

BIAXIAL STRAIN RELATION 78 

Appendix D 

PLATE DEFLECTION PROFILE AND AVERAGE TRANSVERSE STRAIN 80 

Appendix E 

COMPUTER PROGRAM FOR ICE /PLATE INTERACTION 96 

BIBLIOGRAPHY 124 



5 



LIST OF FIGURES 



Number 

1 . 

2 . 



3. 



4. 



5. 



6 . 



7. 



8 . 



9. 



10 . 

11 . 

12 . 

13. 



14. 



Al. 



Title 

Idealized Ice-Load 

Rectangular Plate Dimensions and Load 
Position 

Plate Deflection and Partial Load 
Relationships 

Strain Distribution for Combined Bending 
and Membrane Forces 

Plate Failure Parameters & and 0 for the 
Experimental Plates of References E2 and 
E3. 



Comparison of Plate Deflections to the 
Experimental Results of Loeser (El) 

Comparison of Calculated Plate Deflections 
to the Experimental Results of Day (E2, E3) 

Load versus Absorbed Energy for Day's (E2, 
E3) Experimental Test Plates 

Load versus Absorbed Energy for Loeser ' s 
(El) Experimental Test Plate 

Load versus Absorbed Energy for Arbitrary 
Plates of Tables E12 to El 9 

Plate Failure Parameters 8 and 0 for the 
Experimental Plates of Reference El. 

Plate Failure Parameters 6 and 0 for the 
Arbitrary Plates of Tables E12 to EI 9 

Load/Deflection Relation for the Arbitrary 
Plates of Tables E12 to E14 (Mild Steel 
Material Properties) 

Load /Deflection Relation for the Arbitrary 
Plates of Tables E15 to El 9 (High Strength 
Material Properties) 

Load Coefficient for the Bending Moment 
at the Mid-Point of a Partially Loaded 
Clamped Rectangular Plate (See Equation A3) 



Page 

18 

19 

24 

30 

36 

47 

48 

49 

50 

51 

52 

53 

57 

58 

66 



6 





LIST OF FIGURES (CONTINUED) 




Number 


Title 


Page 


A2. 


Example of the Graphical Solution for 
the Load Coefficient (3 in Equation A3 


67 


A3 , A4 , 
A5 . 


Size of the Load Bearing Surface at the 

Plastic Collapse Load for Values of B/A, 

A/H. and O' /P (See Equation A7) 
o c 


71-' 


Bl. 


Coefficients k and for the Membrane 

Equation for Deflection and the Deri- 
vative of the Equation (See Equations 
B8 and B12) 


77 


Dl. 


Deflection Profile along the Transverse 
Axis of a Partially Loaded Simply Supported 
Plate 


82 


D2 . 


Element of Deflection Profile Curve 


83 


D3 . 


Average Transverse Strain at the Mid-Point 
of the Plate (See Equation D6) 


86 



7 



LIST OF TABLES 



Number 


Title 


Page 


1 . 


Comparison of Calculated Plastic 
Collapse Deflections with Experi- 
mental Results of Reference E4 


42 


2. 


Comparison of Calculated Plastic 
Collapse Deflections with Calculated 
Deflections from Large Deflection 
Theory for a Uniformly Distributed 
Load 


43 


Al. 


Load Coefficients 0 for M of 
Equation (A2) 


65 


A2-A5 . 


Ratio of Yield Stress to Ice Crushing 
Pressure from Equation (A7) 


68-70 


D1-D5. 


Transverse Deflection Profile and 
Average Transverse Strain 


87-91 


D6. 


List of Input Variables for the Computer 
Program for Plate Deflection and Transverse 
Strain 


92 


E1-E19. 


Results of Calculations for Ice/Plate 
Interaction for Various Plates 


99-117 


E20. 


List of Input Variables for the Computer 
Program for Plate Deflection/Load Relation 


118 



8 



NOMENCLATURE 



F > F n’ F o 


-(lbs) 


The force behind the pyramid of ice 


F c 


— (psi) 


Ice -crushing pressure 


q 


"(psi) 


Any uniformly distributed load over a 
previously defined bearing area 




-- (in-lb) 


The bending moment about the longitudinal 
axis 


M o 


-- (in- lb) 


The bending moment at the plastic collapse 
load . 




--(psi) 


The yield stress for a material 




--(psi) 


The ultimate stress of a material 




--(psi) 


Equivalent yield stress of material in an 
elemental strip of plate 




--(psi) 


In-plane membrane stress 




--(psi) 


Total in-plane stress due to the membrane 
force and bending moment 


s,s a ,s i 


--(lb/in) 


In-plane membrane force per unit length of 
plate edge 


u 


-- (in- lb) 

in3 


Energy per unit volume absorbed by the plate 


R 


- - ( in- lb ) 
in^ 


Total restraint energy per unit volume of a 
material in simple tension 


6u 


-- (in/ in) 


The ultimate strain of a material in simple 
tension corresponding to C T u 


6 t 


-- (in/ in) 


The total strain at a- point due to bending 
and stretching 


6 x 


-- (in/ in) 


The strain in the x-direction at a given 
point in the plate 


6y 


--(in/in) 


The strain in the y-direction at a given 
point in the plate 


e z 


-- (in/in) 


The in-plane strain of the outer fibers of 
the plate 
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z,z i’ z i+l 


-(in) 


Partial center deflection as determined 




by membrane formula 


V 


--(in) 


Total plate deflection at any point (x,y) 


W ,w. )W . + i 


--(in) 


Total center deflection for a given load 


y 


(radians) 


Angle between horizontal and midplane of 
plate at a given deflection 


e 


(rad ians) 


Angle of rotation due to bending of a plane 
cross-section through the thickness of the 
plate 


S = W/A 




A non-dimensional failure parameter for 
plate geometry 


0 = U/R 




A non-dimensional failure parameter for 
plate material properties 


k 




Coefficient within the membrane formula that 
is dependent on x and y 


k d 




Coefficient within the derivative of Z that 
is dependent on x and y 
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INTRODUCTION 



Discussion of Problem 



The solutions to two interrelated problems are sought in 
this analysis. The failure of clamped rectangular plates may be 
separated from ice-loading of plates, but the two studies are brought 
together here for the consideration of a special problem in ship's 
structures. The full analysis of hull plating subjected to the various 
loading conditions associated with navigation through areas where ice 
is present brings the two topics together. The first problem is to 
determine the relation between the hull plating and the ice that may 
be forced against the plating. The second problem is encountered as 
the force driving the ice against the ship is increased and the hull 
plating deforms until plate failure or rupture occurs. The total pro- 
cess may be considered continuous, starting with the initial contact 
with the ice, subsequent ice crushing and plate deflection, and finally, 
extreme plastic deformation of the plate and failure. 

The need for this analysis is generated more from the design 
of merchant ships which are to operate in ice filled waters than from 
icebreaker design. Icebreakers are generally very conservatively 
designed because they must be able to withstand ice-loadings of almost 
any magnitude during their lifetime. On the other hand, merchant ships, 
even those with Icebreaker bows, have to balance their design between 
profit and repair costs. Their carrying capacities are greatly affected 
by the addition of thicker steel plates along the ice-belt. Furthermore, 
if there is a chance that the ship could lose part of its cargo through 
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a ruptured plate, the economics of the situation increases. Even more 
important though, is the rising concern over the possibility of oil 
pollution of the arctic seas. In order to prevent oil spills, it has 
been proposed that a tanker be designed with outer buffer tanks filled 
with water. The outer hull would be of thinner plate than the usual 
single hull ship. The separation of the two hulls would be based on the 
maximum deflection that could be expected from the outer plating when it 
failed under ice- loads. 

Idealization of Ice-Loads 

The ice-loads to be considered will be idealized with regards 
to their geometry and point of application. The present design practices 
are based on static loading of the hull plates. The model loads used 
here will be considered continuous but applied in incremental static steps. 
Icebreaker and ice-belt plate scantlings are usually obtained by using a 
uniformly distributed pressure equal to the ice crushing pressure. Since 
the actual ice crushing pressure is not constant and dependent upon 
geographical location, age, salinity and many other factors, it was 
desired to try to cover a wide range of possible loading mechanisms. 

The driving force behind the plate/ice interaction must also be considered. 
It was for these reasons that a pyramid shaped pinnacle of ice was chosen 
as the idealized ice-loading mechanism. 

Type of Analysis 

Plate failure analysis lies between the theory of plasticity 
and a handful of experimental tests that have carried applied loads up 
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to plate rupture. Theoretically, it would be desirable to use the theory 
of plasticity to analyze the plate deflection and failure. This was 
considered, but after reading other works in the field and considering 
the accuracy of the possible results, it was decided to use a more 
empirical approach. The final analysis will be a combination of plastic 
analysis and engineering judgment. This combination of approaches to 
the problem allows a study of the plate/ice interaction within a limited 
time period with sufficient accuracy of the results. 

Load /Deflection Relations 

Most of the experimental work on rectangular plates with large 
deformation has been done under uniformly distributed loads. To obtain 
any data on concentrated loads, it was necessary to go back to a report 
of experiments conducted in 1937. The plate deflections and stresses 
were kept in the elastic range in these tests. The results did provide 
a good starting point for the shapes of the plates during plastic defor- 
mation under concentrated and partial loads. Some of the results and 
observations on circular plates under concentrated loads were borrowed 
for the analysis. 

It was necessary to consider the concentrated or partial loads 
in order to present the complete model of plate/ice interaction. The 
crushing process of ice in contact with the ship’s hull will begin as 
a concentrated load. The ultimate strength of sea ice in compression 
is usually low enough in relation to the design pressures of hull plating 
such that the ice-loads are uniformly distributed before permanent defor- 
mation of the plate begins. However, since the ice crushing process will 
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require some partial loading, large plate deformations may occur under 
smaller, but more concentrated loads. Therefore, partial loading of the 
plate will be analyzed in order to consider all possible extremes for 
the failure analysis of the plate. 

Failure Analysis 

Experimental results of rectangular plate failures are usually 

an added column of data to the experimenter 1 s primary reason for testing. 

The beginning of investigations into large plastic deformation and failure 

is often attributed to Bach. His experiments conducted in 1908 formed 

the basis for Hovgaard's 20% permanent deflection rule.'*' These tests, 

like most of the data obtained for rectangular plate failure, were for 

uniformly distributed loads. Often the comment that the plate failed at 

(T4) 

a given pressure is the only amount of data given. Greenspon, v in 
evaluating experiments by J. W. Day, v 5 attempted to correlate the 
failure data so that plate rupture under a uniformly distributed load 
could be predicted. His recommendation was to establish a limiting 
maximum deflection to span ratio. Using a statistical approach, it could 
be stated that a uniform lateral load giving a deflection to span ratio 
greater than 0.10 could be considered approaching the plate failure 

load. The only other author who discussed plate failure to any extent 

_ (El) 

was Loeser. 

Most of the experiments as reported by the authors listed in 
the bibliography under Plate Deflection-Experimental were helpful in 
establishing a failure criterion. It was decided early in the analysis 
that an energy relation would be used. This enabled correlation of the 

1 See reference T5 , pp. 120-125. 
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failure data that was available. Not all of the data was useful, since 
all material properties must be known in order to properly correlate the 
experiments. Since Day and Loeser included all material properties in 
their reports, the results of their experiments were relied on most 
heavily. 

Outline of Procedure 

The procedure of investigation that follows will continue in 
a similar format as used in the introduction. First, the model of the 
ice-loads will be formed. Then, the plate deflection/load relations 
will be established. At each increment of loading or deformation of the 
plate, the amount of energy absorbed by the plate material should be 
calculated. The relations for this segment of the analysis will depend 
on the load/deflection calculations and will serve to relate them to the 
failure criterion. The final step in the analysis will be the estab- 
lishment of the failure criterion. This will complete the model of the 
ice/plate interaction. 

Ordering of the Bibliography 

It has been previously mentioned that this analysis falls 
between the basic theory of plasticity and empirical relations. Many 
of the assumptions and choices of methods are based on the experiments 
and observations of other authors. Each phase of the analysis can be 
performed by other techniques. However, within the time limits of the 
investigation, all cases were carefully studied so as to make use of the 
best method available. For these reasons, it was felt that the list of 
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references should reflect the degree of influence of each work upon the 
final analysis. Therefore, the references within each subject area are 
listed in approximate descending order of influence or use. 
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PROCEDURE FOR CALCULATING ICE /PLATE INTERACTION 



Initial Assumptions for Ice-Loa d^ 

A four sided pyramid of ice as shown in Figure 1 will be the 
loading mechanism. The manner of loading will be such that it may be 
considered to be static at all times. 

The included angle of the pyramid frustum will not effect the 

load bearing area of the ice on the plate. Also it is small enough to 

cause only the crushing portion of the ice to make contact. That is, 
the sloped sides of the cone will not come in contact with the plate at 
any time. 

The load bearing surface will conform to the curved shape of 

the plate as if the plate in the region of the load is subjected to a 

uniformly distributed load. Thus, load edge stress concentrations will 
be neglected. Also, at no time will the area of load bearing surface be 
less than twice the thickness. Therefore, shear forces will be neglected. 

The load bearing area will be a square of side equal to C. 

If C equals the width of the plate, the load bearing area will be set 
equal to the plate area. Also, at this point, the loading pressure will 
be shifted from a constant ice crushing pressure to a pressure equal to 
the total load divided by the plate area. 

Initial Assumptions for Rectangular Plates 

The plate material will be isotropic and homogeneous. All 
four sides, of the plate will have a clamped edge condition for support 
as shown in Figure 2 and only in-plane stresses are to be considered. 



17 




Figure 1. Idealized Ice-Load 
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Figure 2. Rectangular Plate Dimensions and Load Position 
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Only deflections after the plastic collapse load is reached will be 
considered. A rigid, perfectly plastic material is assumed for the 
collapse mechanism and a biaxial strain condition will be assumed at 
the mid-point of the plate. 

Model of Ice-Loads 

The idealization of an ice- load will be a four sided pyramid 
of ice. The actual applied load will be in the form of a force acting 
on the base of the pyramid. The force will be transmitted to the rec- 
tangular plate in the form of a uniformly distributed load acting over 
an area centered at the mid-point of the plate. The value of the load- 
ing pressure will be the ice crushing pressure except during initial 
loading and after the bearing surface is equal to the plate surface 
area. These exceptions will be discussed later. 

As the driving force increases, the bearing surface of the 
pyramid will increase as the ice crushes. Thus, if the ice crushing 
pressure is considered constant throughout the loading process, the 

bearing surface area will increase directly with the driving force. 

2 

Since the area, C , at any time is equal to the driving force, 
F, divided by the ice crushing pressure, , the length of the side of 
the bearing area will be: 

c ']ff (1) 

This relation will hold after the pyramid first starts to crush until 
the plate is completely covered by the bearing surface. The loading 
prior to the initiation of ice crushing will not be considered, since 
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the plate deflection will be in the elastic range. 

Once the bearing surface equals the area of the plate, the 
applied load, F, will remain constant unless is allowed to increase. 

Thus it will be assumed that ice crushing ceases and will increase 
with any increase in F as: 

P = JL (2) 

AxB 

Plastic Collapse Load 

The collapse mechanism for the plastic analysis is the general 
roof shaped deflection surface. For the transverse cross section a 
three hinge collapse mechanism is used. For any loading, it was assumed 
that the center plastic hinge forms either simultaneously with the two 
edge hinges or at some increased load beyond the formation of the edge 
hinges. 

Since the analysis is for large deflections of plates the 
deflection at the plastic collapse load will be the first to be considered. 
Unfortunately, there is no easy solution for the deflection at this load. 
The problem is further complicated when partial loads are considered. 

Finite displacement relations from the theory of plasticity and a 
membrane solution were considered to be the easiest methods for finding 
the deflections. A similar degree of accuracy may be obtained from 
either method. The membrane solution has been shown to provide good 
representation of the plate deformation for loads well beyond the collapse 
load. Therefore, it was decided to use the membrane solution throughout 
the analysis. 
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As the pyramid of ice is forced against the plate, the pressure 

acting over the initial bearing surface will rise until, the ice crushing 

pressure is reached. During this time, the plate will behave elastically. 

From the time P^ is first reached until plastic collapse of the plate, 

the area of the bearing surface will increase. There will be associated 

2 

with the collapse load F , a constant P , and an area C . In order to 

o c o 

use the membrane solution to find the deflection at collapse, the value 
of must be known. Actually it is easier to use the non-dimensional 
form C/A. 



The moment at the center of the plate was found using an 

elastic solution for clamped plates under partial loads. M = P 

x 

where P is the total load and the load coefficient is given in 

2 

Figure Al for values of B/A and C/A. Then by setting M = M q = C^H /4, 

the yield moment for a rigid, perfectly plastic material, the bearing 

surface area can be related to the yield stress and plate thickness. 

2 

Thus, when P = P^C and both sides of the equation are divided by the 
width of the plate, A, the following non-dimensional equation is 
formed: ^ 



Oo = 4 (3 



'[(SMS)] 



(3) 



From equation (3) the bearing surface area at the collapse 
load may be determined for any plate. Various solutions to the equation 
are given in Figures A3-A5 for given values of B/A and A/H. 



This is equation (A7) of Appendix A. 
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Plate Deflection Calculations 



Once the size is determined for the load bearing surface at 

the collapse load, it is necessary to calculate the deflection of the 

plate at this load. As discussed previously, a membrane solution will 

3 

be used for the deflection calculations. The membrane formula, 

t = 4qA 2 k (4) 

n-3 S 

in general only holds for a uniformly distributed pressure over the 
entire plate. This meant that some geometric relationships were re- 
quired in order to use the membrane formula with partially loaded 
4 

plates. It was assumed that the plate in the region of the partial 
load behaves as a membrane after plastic collapse was reached. Then 
the membrane or in-plane forces along the edge of the loaded area are 
carried through to the edge of the plate. This implies that the slope 
of the deflection profile at the edge of the loaded area is the same 
as the slope at the edge of the plate. Thus the conical shape is 
represented for the concentrated loads and the full membrane shape is 
used when the load is uniformly distributed over the plate. It is 
believed that this provides a good representation of the true deflection. 
The deflections for the concentrated loads may be less than they should 
be, but the accuracy should improve as the load bearing area increases. 

The total deflection calculations consist of three basic 
steps. First, the deflection Z for the loaded portion of the plate is 
found from the membrane formula. For the width of the applied load C 

3 

See Appendix B. 

4 

See Figure 3. 
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Figure 3. Plate Deflection and Partial Load Relationships 



and the pressure P^ , Figure 3 shows the relation of Z to the total 

deflection. The second step is to determine the- slope of the plate 

at the edge of the loaded area. The derivative of the deflection 

profile provides the slope at the point x = (A t C)/2. Finally, the 

total deflection at the mid-point of the plate is found by combining 

the membrane solution with the assumed slope at the plate edge. A 

simple geometric relation will determine the deflection at the edge 

of the loaded area, and the membrane solution may be added to it to 

give the total deflection. The central deflection W is given by: 

W = Z + 1 (A-C) tan , (5) 

2 

where the membrane formula 

Z ■= 4 c 2 k P „ (6) 

7T3 s 

is the deflection of the loaded portion, 

and tan V = d__Z = ? c C kn at x = A-C . (7) 

dx 2S 2 

The constants k and are the factors that depend on x. The values 

for these constants are included in Appendix B for given plate aspect 
5 

ratios . 

Prior to discussion of the procedure for load/deflection 
analysis, the membrane solution for plate deflection must be understood. 
The actual formulation of the deflection relation is given in Appendix 
B. There are two dependent variables in the relation. The membrane 
force, S, which acts along the edge of the loaded area depends on the 
total deflection of the membrane, W, which in turn, depends on S. The 

^ See Figure Bl. 
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usual procedure for obtaining the value of S is to assume a membrane 
deflection and shape. Then a strain relation is used to obtain the 
principal strain in the plane of the plate for the assumed deflection. 

The usual assumed deflection shapes are cylindrical or parabolic. 

In order to model the ice loading condition, the equation for 
the shape of a partially loaded, simply supported plate was used. A 
biaxial strain relat ion which depends on the B/A ratio, the size of the 
bearing surface given by the C/A ratio, as well as the maximum deflection 
of the plate was established for the mid-point of the plate. This made 
it possible to obtain the membrane stress from the stress-strain curve 
for the plate material. Thus, for a given deflection, the membrane 
force per unit length is equal to the membrane stress times the plate 
thickness. 

If the membrane formula (equation 4) is used alone, the applied 
load, q, could be found that would produce an assumed deflection. Since, 
for the model of ice-loads, the loading pressure is assumed constant, 
an iterative procedure must be used to obtain the proper correlation 
between Z and S for a given constant load P . 

For the calculations involving the deflection of a partially 
loaded plate, the total deflection, W, is used instead of Z to determine 
the value of the membrane force, S. This appears to provide a more 
realistic situation, because the strain relation as derived in Appendix 
D depends on both the plate shape and deflection. 
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See Appendix D. 



Load/Deflection Relations 



The procedure for obtaining the deflection, W, for each value 

of the total load, F, will be iterative. The beginning will be based 

on the plastic collapse deflection. After the value of (C/A)^ is found 

from the plastic collapse relation, the membrane solution will be used 

to find the deflection at the plastic collapse load, F^. The first 

approximation for the deflection is found by assuming that the membrane 

stress is equal to the yield stress. Then for the value found for 

from (C/A)^ and the ice crushing pressure, , the first values of Z, 

dZ/dx and W (equations 6, 7 and 5) can be calculated. This is the 

starting point of the iteration for the proper value of for the given 

load F . The successive solutions of W will follow the same format as 
o 

outlined below for the general solution. 

For any given total load F^, there will be a corresponding 

bearing surface area represented by the length of one side C • The 

problem is to find the deflection W. that corresponds to F and C . 
r l - n n 

Let F - F + /\F where the subscript n designates a certain 

no r o 

value of static load, and AF is the incremental value of the total 
load. The subscript o refers to the previous value of applied load. 

For example, will initially represent the collapse load. Then 




will provide the necessary load bearing area to determine the value 

of W . . 
l 

The previous value for membrane force will be used to obtain 
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( 9 ) 



a membrane deflection for C . That is, 



n 



Zj. = k P-c 

nr 3 s 



and it follows that 



and 






W. = Z. 
1 1 



- Pc C n kd 

2S 



( 10 ) 



'1 i (A~C n ) ( dZ \ 

2 \ dx / i 



( 11 ) 



The new values W. and C are used to enter the biaxial strain rela- 
l n 

tions and obtain a new value for the membrane force S^. The next step 

is to return to equations (9), (10), and (11) to calculate Z. 

i+1 

(dz/dx)^-^, and The subscript i-fl is used to distinguish the new 

values of deflection from the membrane force. This notation indicates 
that the deflection calculation is made with the membrane force found 
from the previous deflection, W^. After each new value of W^ is found, 
its corresponding membrane force is determined from the strain 
relations and compared with the previous value of S. Thus, the calcu- 
lation will be concluded when the two values of the membrane force 
compare favorably. 

The above procedure is repeated for each new value of F n . 

When the load is reached which produces a load bearing area equal to 

the area of the plate, that is when C - A, then the loading pressure 

p will cease to be constant and increase with the total load. For 
c 

all subsequent increases in load C/A ~ 1.0 and P c = F n /(A x B). This 
assumption of continuous loading on the plate is just an extension of 
the model of the ice loading beyond the point where the pyramid of ice 
has been fully crushed. 
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Energy Absorbed Relations 



Once the load/def lection relation is found, the energy 
absorbed by the plate should be determined at each value of F in order 
to evaluate the possibility of plate failure. The energy of deflection 
will be the greatest at the edge of the plate except for cases of very 
small load bearing surfaces. For the plate thicknesses usually con- 
sidered for ships 5 the load bearing areas will be large enough to rule 
out the shear effects of the concentrated loads. 

Since the plate failure can be expected at the mid-point of 
the long side of the plate, the energy absorbed at this point should be 
greater than that at any other point. There will be two sources of 
strain energy at the edge of the plate. The initial source will be from 
the bending moment. Then as the membrane forces increase, the portion 
of the total strain due to the bending moment will decrease. The com- 
bination of the strain due to bending and stretching will give the 
total strain at the expected point of failure. 

Using the assumption that plane cross-sections remain plane 
during the bending and stretching, the strain distribution across the 
plate thickness will be as shown in Figure 4. The rotation of the 
cross-section is due to the plate bending, and the location of the 
neutral axis away from the mid-plane is determined by the size of the 
membrane stress. The influence of the membrane stress is equivalent to 
lowering the neutral axis from the geometric center to a new point z 
This shift in axis of rotation of the cross-section depends on the 
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6 Z 




Figure 4. Strain Distribution for Combined Bending and 
Membrane Forces 
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( 12 ) 



ratio of the membrane stress to the equivalent yield stress;^ 

z 1 = h cr g 

2 Op 

where Op = Cf 0 _ = 1-156 O Q for V = 0.5 (13) 

Vl- l/ - 12 2 

This relation will hold true for values of (J greater than 0' since 

H P 

at that point the neutral axis will lie "outside" of the material and 

all strains across the thickness will be in tension. 

The strain of the outermost fiber at the mid-point of the 

long side will be represented by 6 . From Figure 4 

2 1 5 

&Z = e s + | tan 0 (14) 



Since 6 

s 



z 1 tan 8 = H (T s tan 8 
2 On 



e 

2 



= H tan 8 

2 



i) 

a p J 



(15) 



Most experiments have, shown that the longitudinal strain, 

6 at the mid-point of the long side of the plate is very small 

compared to the transverse strain at that point. Thus £ was taken 

5 2 

as the principal or transverse strain at the point of expected failure. 
Again using a biaxial strain relation as given f n Appendix C, the 
uniaxial strain is, 



6 



t 



6 

1- 



(16) 



where V is equal to Poisson's Ratio. 



Young in reference E5 recommended the use of Qf i n the equation 
for the shift of the neutral axis due to the addition of membrane 
stresses . 
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In order to find the stress at the expected point of failure, 
the uniaxial strain was used in the constituat ive relation for the 
plate material. The total stress will be designated 0 and may be 
found for the strain 6 t from the stress-strain diagram of the material 
in simple tension. 

The energy absorbed per unit volume at the mid-point of the 
long side will then be represented by: 

u = e t cr t lb s J (17) 

This value may be found at any point in the loading sequence if it is 

assumed that 6, the angle of rotation of the cross-section, is equal 
to y the angle of the slope of the deflection profile at the edge 
of the plate as shown in Figure 3. This is a reasonable, as well as 
a convenient assumption, since the slope at the edge is determined 
during the deflection calculations. Thus, for each value of F and its 
corresponding W, 

tan 0 - dZ (18) 

dx 

Then 6 t may be rewritten as 

6. = H / dZ \l CTs + l\ (19) 

2(1-V 2 ) [ dx j\ CXp ) 

Also, the stress-strain relations for the plate material that are 
required for the deflection relations may be used to find (J ^ . There- 
fore, for any loading and corresponding deflection, the energy absorbed 
by the plate at the mid-point of the long side may be determined. A 
comparison of this energy to the restraint energy of the material will 
lead to a failure criterion. 



32 



Failure Criterion 



When a plate is drawn in biaxial tension, the strain rate 

becomes unstable when the principal stress reaches the ultimate strength 

of the material. The plate begins to decrease in thickness rapidly like 

(Pi) 

a tension specimen begins to neck at the ultimate stress. For this 

reason it was decided to use the ultimate stress of the plate material 
to establish the failure criterion. The ultimate stress will be defined 
as the stress at which the slope of an engineering stress-strain diagram 
is equal zero. The area under the stress-strain curve for all strains 
less than the strain at the ultimate stress will represent the stability 
of the plate material. 

A restraint energy relation will be defined as 

( 20 ) 



r = e o 

u u 



in- lbs 
~Tn3 



where O is the ultimate stress and £ is the ultimate strain, 
u u 

Both values are defined at the point of zero slope of the stress-strain 



curve . 



The term restraint energy is used since R represents the 
amount of energy that might be absorbed by the material before insta- 
bility begins. A more accurate measure of the ability of the material 
to absorb energy would be the area under the stress-strain curve. Since 
this area is usually difficult to obtain readily, the R relation is used. 
It provides a reasonable representation of the area under most stress- 
strain curves. Furthermore, it is more desirable to obtain a relation 
that will give a consistent measure of the restraint energy rather than 
an accurate one. 
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The two parameters that best correlate the minimum amount 
of failure data that is available are the deflection to span ratio, 

W/A, and the absorbed energy to restraint energy relation U/R. These 
non-dimensional parameters represent virtually all of the factors 
involved in the plate problem. Several other parameters were considered, 
but there did not appear to be any relation between them and the failure 
data. The W/A ratio helps to bring together all the geometric proper- 
ties of the plate. The material properties are well represented by the 
U/R ratio. 

Figure 5 is a plot of the two non-dimensional parameters for 
most of the rectangular plate failure data that is available. The curves 
are drawn through data points obtained by the calculation procedure 
described in the previous sections. The failure points for each plate 
are designated in two positions. The calculations for the deflections 
of the plate did not correspond directly with the experimental data as 
given in references E2 and E3. In general, the calculated deflections 
were slightly lower than the experimental deflections for a given total 
load. The possible reasons for this discrepancy will be discussed later. 
It was felt that both failure deflection and failure load should be 
indicated in Figure 5. Thus, the first set of failure points were found 
by using the experimental deflection at the mid-point of the plate to 
determine the W/A ratio at failure. The second set of points were ob- 
tained by using the failure load from the experimental data. For each 
rupture pressure, the total loading force F was calculated and from 
Figure 8, the corresponding U was found. The U/R values for the failure 
loads are shown on Figure 5 as the lox^er set of failure points on the 
curves. 
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The lower set of failure points was chosen to represent the 
failure data. Although the lower set provides the more conservative 
choice for a failure criterion, the major reason for using this set is 
that the failure data corresponds better to the deflection and absorbed 
energy calculations described in the previous sections. 

A line can be drawn on Figure 5 connecting most of the fail- 
ure points. With the exception of the plate designated IB this line 
represents a minimum relation of W/A to U/R for plate failure. If it 
is moved slightly to the left, the failure point for plate lB will also 
lie above the line. This new position of the limit line will be used 
as the failure criterion. 

If 6 - W/A and 0 = U/R, then failure may be expected beyond 
the line defined by 

6 - .278 - 0.182 0 (21) 

Greenspon in reference T4 used essentially the same failure 
data and stated that failure should be expected for loads beyond the 
pressure that produces W/A =» 0.10. This limit as shown in Figure 5 is 
quite conservative for 0 <0.9, but for plates with 5 <0.10 and 0 > 0.9 
the limit may be too large. 

The energy relation was chosen for this analysis in order to 
provide a parameter that will accurately represent all types of loading 
conditions. Unfortunately, the available failure data is for loads 
uniformly distributed over the entire plate. In order to properly 
relate the failure for partial loads and uniform loads, equivalent 
absorbed energy relations are used. Thus the failure criterion should 
hold for all but concentrated loads with small bearing surfaces. 
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Summary of Meth o d of Analysis 



The method of calculation of the deflection and the absorbed 
energy for any plate has been incorporated into a short computer program 
The description of the program itself is in Appendix E. The only dif- 
ference between the procedure summarized below and the computer program 
is the method of determining a stress for a given strain. A generalized 
equation for the plastic portion of the stress-strain curve is used in 
the program. With the proper choice of variables the generalized stress 
strain equation can represent the actual stress-strain curve of a given 
material with reasonable accuracy. There are several variations to the 
general equation that may be used. The relation described in Appendix 
E appeared to be one of the easiest to use. 

The computer provided the speed and accuracy to carry out 
the calculations for many plates. The hand calculations are not diffi- 
cult, but they are quite tedious. With the use of the tables and curves 
presented in the Appendixes, the time for calculation can be shortened 
considerably. The step-by-step calculations for a plate of known geo- 
metric and material properties and a given ice-crushing pressure are; 

1. For the values of B/A, A/H 5 and G /P enter the curves 

o c 

given in Figures A3, A4 or A5 to obtain the value of C/A at the collapse 

2 

load. Then calculate F - C P . 

o c 

2. Calculate Z = 4 C^k P c , and dZ = T r C krf , 

IX 3 S dx 2S 

and W = Z -1 JL (A-C) _d_Z where k and k^ may be 
2 dx 

obtained from Figure Bl for B/A. Also the membrane force for the first 
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approximation will be assumed to equal the yield stress times the 
thickness . 



3. The quantities W/A, B/A, and C/A are required to find 

the value of the average transverse strain at the mid-point of the plate. 
Tables of values of the average transverse strain and representative 
curves are given in Appendix D. 

4. The new value of membrane stress will be found from 
the stress-strain curve for the plate material for the average trans- 
verse strain found in step number 3. Then = 0 # s H. 

5. Steps 2 5 3, and 4 are now repeated as described on page 
28 until the new value of is within an accepted difference from the 
previous value. (An absolute difference of 25 psi is used as a tolerance 
in the computer program) . 

6. For the final values of CX S and dZ/dx from step 5, cal- 
culate 6 - H ( dZ Cfg + ^ J . 



2 ( 1 - V 2 ) 



dx 




Or 



7. Again enter the stress-strain curve for the plate material 
to find CJ t which corresponds to 6 t . 

8. Calculate the energy absorbed at the point of expected 

failure. U = 6 t CT t 

9. Use the value for 6 u as defined on page 33, and calcu- 
late R = 6 u Cf u . 

10. Using & - W/A and 0 = U/R 5 solve the failure criterion 
equation FC = $ + .182 0 - .278. If FC >: 0 5 then failure can be 
expected for any additional loading. 
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11. If FC < 0, increase the total load F by the load incre- 
ment A F which may be any predetermined amount. Obtain a new value 
for C =•• Vf/P c 

12. Repeat steps 2 through 11 to find the deflection and 
absorbed energy values for the new loading condition. 

It is not necessary to evaluate the possibility of failure 
at each load increment. After one or two times through the procedure 
the approximate point of failure may be apparent and the solution 
method can be shortened. 
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RESULTS OF APPLICATIONS OF THE METHOD OF ANALYSIS 



Plastic Collapse Load and Deflection 



The method for determining the total load for plastic collapse 
depends on a given ice crushing pressure. The area over which this 
pressure must act to cause collapse is determined as outlined on page 22. 
If this method is reversed so that for a given C/A and A/H a certain 
<V P C is found, then the required pressure for collapse will be known. 
When this pressure is compared to elastic experimental data for plates 
of the same size and material properties, the collapse pressure should 
be considerably greater than the elastic loading pressure. It is very 
difficult to compare deflections, although the deflection should be 
increased by an amount similar to the loading. 

The plastic collapse load and deflection were calculated for 
plates of the same size and material properties as those used by Sturm 
and Moore . The size of the bearing surface was fixed at 2 M x 2 M 

and 12 M x 12". Then the plastic collapse load was found from Figures 
A3-A5. In this case P c , which has been called the ice crushing pres- 
sure, is the collapse load. Finally, using these values for C and P c 
the collapse deflection, W G , was determined by steps 2-5 of the previ- 
ous section. The results of the calculations for three of the plates 
are shown in Table 1 along with experimental data from reference E4. 

A second comparison of the plastic collapse loads and deflec- 
tions was made with the results from the large deflection theory of 
plates. This provided a means of judging the accuracy of the membrane 
solution method for calculating the plastic collapse deflections for 
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uniformly loaded plates. 

Reference E8 compares many of the solutions to the large 
deflections of plates to experimental results. Most of the relations 
are given in the form of curves for the non-dimensional deflection, 
stress, and pressure ratios. The deflections of three plates were 
calculated from the corresponding pressure ratio. The plates were of 
the same size and material as those given in Table 1. The collapse 
pressure obtained from Figures A3-A5 for C/A =1.0 was used in the 
pressure ratio, 

P. R . = Po aA (22) 

16EhA 

For example, Figure A3 shows that (T 0 /P c = 215 for an aluminum 
(E = I0xl0 6 psi) plate I"x48"x48" with 0 = 38800 psi and C/A = 1.0. 

Thus P = p G = 180.5 psi. Using the value in the pressure ratio, a 
deflection ratio W/lI =0.9 was obtained from Figure 24 of reference E8. 
The same procedure was followed for the other two plates. The results 
from the large deflection theory and the membrane solution for collapse 
deflections and loads are shown in Table 2. 



The Membrane Solution for Plate Deflections 



It was relatively easy to determine the accuracy of the mem- 



brane solutions for the plate deflection. The plate dimensions and 

mate 
(El) 



(E2 E3) 

material properties of the test plates used by Day 5 and Loeser 



were put into the computer program to determine the load/deflection 
relationships. The results of the calculations as well as the input 
data are shown in Tables El to Ell of Appendix E. It should be noted 
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" LOAD BEARING SURFACE 
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Table 1. Comparison of Calculated Plastic Collapse Deflections with Experimental 
Results of Reference E4. 
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Table 2. Comparison of Calculated Plastic Collapse Deflections with Calculated 

Deflections from Large Deflection Theory for a Uniformly Distributed Load. 



that the input variable for ice crushing pressure for these calculations 
is only the loading pressure that is used for the first . calculation. 

The pressure used for any additional loading is equal to the applied 
force divided by the plate area. It was necessary to follow this pro- 
cedure because all of the test data is for uniformly distributed loads. 
Actually the solution method is the same as if an ice-load is used and 
C/A = 1.0. In this way the computer program may be used for any type 
of lateral load except highly concentrated loads. 

Also, the first calculation for deflection given in Tables 
El to Ell does not represent the collapse deflection. The starting 
point used for these calculations was chosen for convenience in com- 
paring the results to the respective experimental data. For example, 
the rupture pressure given in Day f s experiments was divided by five and 
the result was used as the input for ice crushing pressure in the com- 
puter program. Then the increment of applied force & F was set equal 
to the ice crushing pressure times the plate area. 

The results of the calculations are shown in Figures 6 and 7. 
The experimental results of Loeser and Day are also plotted in each 
figure . 

Energy Absorbed Relation 

It was of interest to plot curves of total load F to the 
energy absorbed U for three sets of data. The F versus U relations 
for the calculations on Day's test plates as given in Tables El to E7 
are shown in Figure 8. Similar curves for calculations on Loeser's 
test plates (Tables E8 to Ell) are plotted in Figure 9. The results 
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of calculations made on arbitrary plates under idealised ice-loads 
are given in Tables E12 to E19, and the F versus U curves are shown 
in Figure 10. 

The calculations for Loeser f s and Day 1 s test plates have 
been previously described. The results shown in Tables E12 to E19 
are for four plates with the same B/A ratio, but of varying thickness. 
Also, two sets of material properties were used for each plate. The 
calculations follow the format as described in the method of analysis 
and Appendix E. 

A particular case is encountered in the calculations for 
the arbitrary plates that was bypassed for the calculations of Loeser's 
and Day f s test plates. The assumption that the ice crushing pressure 
changes when C/A - 1.0 has been discussed previously. The analysis of 
the arbitrary plates required this change. This can be seen by tracing 
down the C/A column in Tables E12 to E19. For all but the thicker 
plates C/A starts at the value for (C/A) c and increases until C = A. 

For square plates the ice crushing pressure will change to the total 
load divided by the plate area. But for plates where B/A > 1.0, this 
change can cause a decrease in the crushing pressure rather than the 
desired increase. As a consequence the deflection will be smaller than 
the previous deflection. In order to help eliminate this unnatural 
discontinuity, it was assumed that the loading surface continued to 
grow after C/A = 1.0. That is, the loading pressure was kept constant 
until not just the width of the plate but the entire plate was covered. 
Then for any further increases in the load, the pressure increased as 
the force divided by the area of the plate. This scheme would be the 
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natural procedure to follow if the calculations were being done by 
hand but it required a few extra statements in the computer program. 

Application of Failure Criterion 

The relation between $ - W/A and 0 = U/R has been plotted 

in Figure 5 for the results shown in Tables El to E7. The failure 
criterion was established from these curves and the failure points 

/DO 

given by Day 5 . The establishment of the 6 versus 0 curves is 

discussed on pages 33-35. 

For the W and U values given in Tables E8 to Ell, the corres- 
ponding 8 f s and 0’s were calculated. The relation between the two 
failure parameters is shown in Figure 11. These curves represent the 
experimental work of Loeser 

Similarly, Figure 12 is a plot of S and 0 for the calcula- 
tions made on the arbitrary plates given in Tables E12 to E19. The 
arbitrary plates by their sizes and properties represent typical hull 
plating on ice-strengthened ships. Thus, the curves of Figure 12 
give a good indication of the possibilities of failure of these plates 
under the idealized ice-loads. 
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Figure 6. Comparison of Plate Deflections to the Experimental Results of Loeser (El) 
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Figure 7. Comparison of Calculated Plate Deflections to the Experimental Results of Day (E2,E3) 




Figure 8. Load versus Absorbed Energy for Day's 
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Figure 12. Plate Failure Parameters 8 and 0 for the Arbitrary Plates of Tables E12 to E19 



DISCUSSION OF RESULTS 



Load/Def lection Relations 

Figures 5 to 11 and Tables 1 and 2 are presented so that a 
comparison of the various portions of the method of analysis may be 
made with associated experimental data. The degree of accuracy of each 
step of the analysis is shown in the tables and figures. 

Although the data presented in reference E4 was published in 
1937, it proved to be the most useful experimental data on rectangular 
plates under partial loads. The comparison between Sturm and Moore 1 s 
experiments and the calculated values for collapse deflection and load 
given in Table 1 is very crude. The experimental data is for the 
elastic range only, and the calculated values assume that the plastic 
collapse load has been reached. However, the experimental, data does 
provide at least a lower bound to the calculated values. 

There are two general points that should be recognized from 
Table 1. The relation between the deflections and the B/A ratio that 
is indicated by the experimental data holds for the calculated values. 
That is, the experimental deflections increase as B/A increases and 
the same holds true for the calculated values. A similar comparison 
between deflections and plate thicknesses shows that the experimental 
deflections increase as the plate thicknesses decrease. However, the 
calculated deflections do not follow the same trend as the experimental 
results. It seems that the thinner plate should have the greatest 
deflection as it does in the experimental data. Although only relative 
comparisons have been made, the failure of the membrane solutions for 
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the collapse deflections to compare favorably with the experimental 
results should be considered a limitation of the method of analysis. 

A comparison of the plastic collapse deflection obtained from 
the membrane solution with the deflections calculated from large deflec- 
tion theory may be made from Table 2. In general, the membrane solution 
deflections are lower than those obtained with the same pressure from 
large deflection theory. For this case of uniformly distributed loads, 
both sets of calculations show that the deflection increases as the B/A 
decreases. Also, the deflection for the large deflection theory remains 
constant for the three plates considered. Again the relation between 
deflections and plate thicknesses does not compare favorably. The 
membrane solution maintains approximately the same W/H ratio at collapse 
for all thicknesses. If the results from the large deflection theory 
are used to define the plastic collapse condition, then the limitation 
in the method of analysis as discussed for Table 1 is again encountered. 

The accuracy of the membrane solution may be easily seen in 
Figures 6 and 7. The calculated deflections are either completely above 
the experimental deflections or completely below. In all cases, if the 
plate has strength properties (see Tables El to 11) greater than mild 
steel ( 0' o = 38000 psi) the deflections lie below the experimental 
data. For the higher strength plates, the deflections are all below 
the experimental results. The one exception is Day’s test plate 2A in 
Figure 7. It has a yield strength slightly above the mild steel plates 
and the curves for the calculated and experimental deflections coincide. 

Since the shapes of the calculated and experimental deflection 
curves agree very well, the only major discrepancy is the difference in 
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the deflections as noted above. The reason for the differences is the 



solution method for obtaining the membrane force for the membrane for- 
mula. This same reason may be cited as the cause for giving lower 
values of deflection for high strength materials and higher deflections 
for lower strength materials. The generalized stress-strain curve used 
in the computer program usually gives values of stress that are lower 
than the actual value for the same strain. This will in turn produce 

larger membrane deflections since the membrane force is S = (X 0 H. 

s 

Furthermore, for two plates of approximately the same thickness, the 
membrane stress will be greater for a plate made of high strength metal 
than one with a leaver yield strength. Consequently, the membrane de- 
flections for the high strength material will be less than those of 
materials with a lower yield stress. 

For the general discussion of the application of the method 
of analysis, Figures 13 and 14 illustrate what the load/deflection 
relationship is for hull plating under the idealized ice-loads. The 
curves and the respective results of the computer calculations given 
in Tables E12 to E19 show that for small load bearing surfaces the 
collapse deflections appear to be excessive. The collapse deflections 
for the large bearing surfaces are usually smaller than what might be 
expected. For the uniformly loaded plates, the collapse deflection is 
not necessary for the load deflection analysis. Thus, for plates that 
require a collapse load greater than that given by step 1 of the method 
of analysis, the first value for deflection given in Tables E15, E18, 
and E19 does not represent the collapse deflection. 

For the partial loading conditions there is no direct relation 
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Figure 13. Load/Deflection Relation for the Arbitrary Plates of 
Tables E12 to E14 (Mixd Steel Material Properties) 
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Figure 14. Load/Def lection Relation for the Arbitrai-y Plates of 
Tables E15 to E19 (Hign Strength Material Properties) 
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between the deflection at one value of the applied force and the next. 
Thus, once the uniform load is reached, the deflection should be within 
the same accuracy as the plates shown in Figures 6 and 7. This fact 
alone should help to substantiate the loading process. In all cases, 
the deflect ion/ load relation, in general, follows a continuous curve. 
There are exceptions to this, but they are caused by the change over 
from the constant ice crushing pressure to the uniformly distributed 
load. If, for example, the deflection for the first partial load is 
greater than the deflections for subsequent loadings, the deflection 
will be decreasing as the total load increases. As can be seen in 
Figures 13 and 14, the deflection is never less than the previous value. 
Therefore, if the accuracy of the calculated deflections for the uni- 
formly distributed loads can be accepted, then the calculated deflec- 
tions for the partial loadings should be of reasonable values for use 
in the failure analysis of the plates. 

Absorbed Energy Relation and Failure Criterion 

There is one basic purpose for calculating the energy absorbed 
by the plate during the loading operation. The calculated energy pro- 
vides a means of determining the possibility of plate failure on any 
increase in the load. There are several methods of determining the 
amount of energy absorbed by a plate. The method presented provides a 
measure at the point of expected failure. 

It has been found in experimental tests of dynamically loaded 
clamped plates (see references E19 and E20) that a linear relation exists 
between the measure of appl i^d force and the measure of the energy 
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absorbed by the plate. The relation between F and U in Figures 8 and 
9 is clearly linear for all plates. The results shown in Figure 10 for 
the arbitrary plates under the idealized ice-loads are not linear. The 
nonlinearity is due to the changing size of the bearing surface whereas 
the curves of Figures 8 and 9 are for a uniformly distributed load 
throughout the loading process. 

It is interesting to note that the only difference between 
the two curves in Figure 10 is the material properties of the plates. 

The F versus U curve for each of the four plate thicknesses considered 
coincides with the other plates of the same material properties. This 
same general relation appears to hold for the plates shown in Figures 
8 and 9, although B/A ratios and plate sizes differ in these figures. 

Thus, based on the results of Figures 8 and 9, and the consistency of 
the results shown in Figure 10, the absorbed energy relation, U, should 
provide a good measure of the possible failure of a plate. 

The calculated values for 8 and 0 from the data given in 
Tables E8 to Ell are plotted in Figure 11. The values of U representing 
the maximum load that was applied to each of the experimental plates by 
Loeser are also indicated. With the exception of plate number 

five, all of the maximum load points lie on or above the failure cri- 
terion line. Plate number five failed due to a weld defect, but the 
other plates did not fail in Loeser f s experiments. Thus, there appears 
to be at least a rough correlation of the failure criterion for two 
different experimental results. 

Figure 12 shows the final application of the method of analysis. 
If the line of the criterion is accepted for the smaller values of 8 
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then the thicker plates can be expected to fail shortly after their 
collapse deflection is reached. This is actually what would be 
expected for the very stiff plates. However, it seems that the line 
of the failure criterion should include larger values of 0 for the 
smaller values of S . This is just a speculation, and the presently 
defined failure criterion appears to provide a reasonable estimate of 
the failure load of a plate. 

CONCLUSIONS 

The method of analysis for large deflections of a clamped 
rectangular plate under idealized ice-loads provides a reasonable means 
of predicting the behavior of a plate. Furthermore, the solution 
method provides a simultaneous check on the possibility of plate fail- 
ure. The calculations are tedious, but not difficult, and the use of 
the computer program will greatly aid any extensive failure analysis. 

A failure criterion was sought and obtained. Furthermore, 
the criterion is not limited to plates under ice-loads. It should be 
equally applicable to the analysis of a plate under any lateral load 
except in cases where shear effects are large. 

RECOMMENDATIONS 

The major recommendation is quite obvious, that is, to conduct 
tests to provide more data for the failure criterion. Other areas that 
should be investigated are other possible solution methods for the 
load/deflection relations and the energy absorbed relation. These areas 
should be included in any failure tests since they are necessary for 
the establishment of the failure criterion. 
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Appendix A 



PLASTIC COLLAPSE LOAD 



The elastic moment at the center of a partially loaded simply 
supported rectangular plate is given by Timoshenko in reference Tl. The 
moment about the longitudinal axis at the mid -point is, 

M 0 P (Al) 

X 

where P = qC^ and (3 is given in Tables 20-22 of reference Tl. The 
values of 0 depend on C/A and B/A, and those values given in Tables 
20-22 decrease in accuracy as C/A decreases. Values for (3 for the 
smaller values of C/A can be found from the equation for M given on 



(A2) 



p. 160 of reference Tl. That equation is: 

M = P T (2 log 4A + X - 1.571) (1+V ) + /JL { 1- V) 

x 8n L nc n 

The parameters X and fX depend on B/A and the values for them are 
given in Table 27 on p. 161 of reference Tl. The values of (3 calcu- 
lated from equation (A2) for C/A less than 0.20 are given in Table Al. 

In order to use equation (Al) for a plate with a clamped edge 
condition, the coefficient (3 must be adjusted, since the values given 
in Tables 20-22 of reference Tl and Table Al are for a simply supported 
plate. The correction term (3 \ may be found on p. 206 of reference Tl 
for various B/A ratios. Thus, the maximum moment at the center of a 
clamped rectangular plate under a partial load is, 

<M ) max = ( p + Pi) P = §'p (A3) 

The factor is plotted against C/A for three values of B/A 

in Figure Al. Two steps were necessary in finding the proper value 
of (3 7 . First, (3 was plotted against C/A for the values of 0 given 
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in Tables 20-22 of reference Tl and for those values of (3 given in 
Table Al. Then for each B/A, a single curve was formed which combined 
the two solutions of 0 .An example of the graphical solution 
for (3 is shown in Figure A2 . The second step required adding the 
constant (3 ^ to the new values of (3 obtained graphically to complete 
the solution of . 

Once the elastic solution for the moment at the mid-point of 
a clamped rectangular plate was established, the plastic collapse load 
was found by simple plastic analysis. The plastic moment at a cross- 
section through the thickness of a rigid, perfectly plastic plate is 

M 0 = Co H 2 (A4) 

4 



Setting (M ) max = M 0 the plastic collapse load becomes, 
x 



Po = Oo H 2 
4 p' 

2 2 

Since P G = q 0 C - P c C for an ice crushing pressure equal to P c 
c 2 = (To H 2 

p c 



2 

By dividing through by the. width of the plate, A, equation (A6) may 
be non-dimensionalized. 




(A5) 



(A6) 



or 



C = _H_ / CTo ' 
A 2A y 0' P c 





(A7) 



The desired relation between the load bearing surface and the ratio of 
material and ice properties for given plate geometry may be obtained from 
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equation (A7). 

The parameter (X 0 /P c for values of B/A, C/A, and A/H has been 
tabulated in Tables A2, A3, and A4. From these values the curves of 
A/H versus Cf Q /P c for various values of C/A have been drawn in Figures 
A3-A5. It was found that this was the best way to present the data in 
order to make use of them in the plate analysis. For a rectangular 
plate of given dimensions and yield stress being loaded by ice with a 
given crushing pressure, (C/A) Q , the size of the load bearing surface 
at the plastic collapse load may be found for A/H and CT 0 /P c from the 
curves of Figures A3-A5. 
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Appendix B 



DERIVATION OF MEMBRANE FORMULA FOR PLATE DEFLECTION 



The solution for the deflection of a uniformly loaded mem- 
brane is given in several references. (See for example, references 

Tl, T4, and El). The method presented here is the one given by 
(El ) 

Loeser. Since the method has been developed elsewhere, only the 

major points of the solution are given. 

The differential equation of a uniformly loaded rectangular 
membrane is: 

J> 2 z . a2z , - -q ( B1 ) 

Vx 2 9 y 2 s 

where q is the lateral loading pressure and S is the membrane force. 

For the plate dimensions given in Figure 2, the boundary conditions are: 



Z - 0 at X " 0, X=A, y =+ B 

2 

The assumed solution for equation (Bl) is the series: 



oo 



Z * 



b sin nftx Y 
n — : n 



(B2) 



n«l,3,5 

The t> n 's are constant coefficients, and Y n 's are dependent upon y only. 

After differentiating twice with respect to x and y and substituting 

for Q Z Z and 3 Z Z , equation (Bl) becomes: 

C>y 2 ~ 



CO 

T 



b n sin n ftx 



n=l ,3,5 



A 



Yn + Y n J = - f 



(B3) 
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Expansion of the right hand side of equation (B3) into a 
Fourier series and using constant value for q between 0 5 x < A 
results in the following expression: 



oo 



£ = (- 1 )' 

« ns z — ' „ , > 

n = 1,3,5 



n-1 



sin n TTx 



(P4) 



Setting equations (B3) and (B4) equal to each other gives 

n-1 



< fl V 



A* 



n 



(-iy 



(B5) 



s nrrb 



n 



Solution of this ordinary differential equation for the conditions 
that Y n is symmetric about y = 0 and that = 0 when y = 4- B/2 yields. 



n-1 

. ia A . 2 _{ri ) 2 



3 3 

S n n b 



cosh nTT B _ cosh n rty 
cosh n TX B ( 2A A 

2A 



(B6) 



Then the membrane deflection may be obtained by substituting equation 

(B6) into equation (B2). 

oo 

cosh n/T y \ 

A j sin n TTx (B 7 ) 

cosh n TZ B J A 

2A / 

If the constant terms are removed and the series is defined as k, the 
membrane deflection equation is: 



Z = 



n-1 



n 



= 1,3,5 



4qA (-1 )' 

3 3 

s n rr 



and 



Z 



4 q A 2 k 
S n 3 



OO 




n-1 




sin n ft x 
A 




cosh n ft y \ 

a_ 

cosh n ft B / 
2A~ / 



(B8) 



(B9) 
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The values for k have been plotted in Figure Bl for three 
values of x/A. The use of I hose curves eases the burden of computation 
for the deflection. For example, the value of the central deflection, 
W, of a uniformly loaded plate will require the value of k correspond- 
ing to the B/A ratio of the plate and x/A » 0.50. This same value of 
k will hold for all values of; q and S. 

The slope of the membrane surface is obtained by differen- 
tiating equation (B8) with inspect to x. Thus, the slope at the mid- 
point of the long side of the plate is: 



oo 



dZ 4qA 



n-1 



9* srt 2 rn /3 

* J i 



lz±ll 



n 



Then for y = 0 and x = 0 or * = A, 

0-0 

dZ _ 4qA^ 

c> x srr 2 

OO 




1- 



n - 1,3,5 



cosh n TX y 

A_ | cos nltx 

cosh n 7T B 
2A 



(BIO) 



1 / - 1 
n^ ( cosh n fr B 
2A 



Noting that ~ j* 7 o r J ~ , equation (Bll) becomes 

n 55 8 



n 



(Bll) 



oo 

9Z qA_ 4qA \ 

3* " 2S ' sn „4-T 



n = 1,3,5 " 



2 cosh nHB 



2A 



or 



0 Z _ _qA k j at x =*» y = 0 

Sx 2S 



(B12) 



where 



k d = i- 8 



TX2 



OO 



n = 1 •> 5 n 2 cosh H-B 1 

,J>:> 2A 



The coefficient is also y, i ve.f in Figure Bl for values of B/A. 
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A ppend ix t; 

BIAXIAL STRAIN INFLATION 



Assume that a given element ol the plate has only normal 
acting on each face. 




forces 



For the case of uniaxial loading let S = S =0 Then the 

y 2 

strain components along each axis will bn 

6 xl* e l’ ®yl " ve xl* c z i = _y e xl _(cl) 

where 6^ represents the principal strain j_ n ^he x direction. The 

strain components in the y and z directi»mg are obtained by considering 

S = S = 0 and S = S =0 respectively, Thi\<* 
x z x y J inus, 



II 

CN 

>> 

(i) 


G 2 5 


^x2 


= 


C z2 = -V6 y2 


(C2) 


and € z3 = 


6 3’ 


G o 

y3 


G z3’ 


C x3 “ 6 z3 


(C3) 


If all 


the 


forces 


are acting, 


1 he total strain in the 


x direc- 


ion will be, 

= 


6x1 


+ ^ x2 


+ G x3 = 6 


1 - V£ y2 ‘ V e z 3 




e x = 


e l 


- V ( 


6 y2 + 6 z3 } 




(C4) 
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For plates it is generally assumed that G is sraall compared to the 

other strains. The total strain, in the x direction for £ =0 will be 

z 

6 X = 6l “ V' £ y2 (C5) 

Similarly, the total strain in i [ ie y direction for £ =0 will be, 

z 

6 y = £ 2 < C6 > 

For the case where the in-plane strain is uniaxial in the x direction 

€ ~ 0 and , hence , 6 = ^ € . = V 6 . From equation (C2) G = £ 

y 2 x i i y 2 2 

= l/6 1 - Substituting into equation (C5) 6^ = ( 1- V 2 ) 6 L . (C7) 

This relation will be useful in obtaining the stress at a given point 
where the. strain in the x direction is known, and the strain in the y 
direction is zero. The stress ran be found from a stress-strain diagram 
for simple tension by using 

€ i " -- x 

i- V 2 - (C8) 

A more general relation is required when there is strain in 
both the x and y directions. ll £ r A 0, the strain in the x direction 

y ? 

may be found by combining equations (C5) and (C6) . 

= e l _VG y2 = G l " V6 2 

= - C ( £ y - Cep 

6 X = c 1 " V 2 ) e 1 -v e y (C9) 

When used in the form 

6 -I ~ £ x -f- V G y 

i -v 2 (do) 

the stress in the x direction mny be found from a stress-strain diagram 
for simple tension. 
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Appendix D 



PLATE DEFLECTION PROFILE AND AVERAGE TRANSVERSE STRAIN 

In order to use the membrane formula for the calculation of 
plate deflections, a relation between the deflection and the membrane 
force must be established. This is accomplished by assuming a plate 
deflection profile. Then, for a given deflection, the average strain 
within the plate can be calculated. From the stress -strain curve for 
the average strain, a value for the in-plane or membrane stress can be 
found which is converted to the membrane force by multiplying by the 
plate thickness. 

For this analysis, it was desired to use an assumed deflec- 
tion profile that would represent the various loading conditions. 

Thus, the solution for a partially loaded simply supported rectangular 

(Tl) 

plate as derived by Timoshenko was selected. The general expres- 

sion from page 138 of reference Tl for the deflection at any point 
along the x axis (see Figure 2) is: 




cosh(0C m -2 7r m ) 



+ 7T sinh( <X -2 V ) + OC 

m mm m 




sin m ftx (Dl) 

A 




Since only the shape defined by this expression is desired, the factors 
in front of the summation sign were dropped. Then, in order to obtain 



a non-dimensional profile, each value of deflection is divided by the 
deflection at x = A/2, the mid-point of the plate. The result is a 
series of shape coefficients which depend on B/A, C/A, and x/A. 

The calculations for w /w for various values of B/A, C/A, and 
x/A were made much easier with the aid of the digital computer. The 
short program used to make the calculations is included at the end of 
this appendix. Table D6 is a list of the input variables required 
for the program. This same program is incorporated into the overall 
ice/plate interaction program discussed in Appendix E. The output 
from the plate deflection profile and transverse strain computer pro- 
gram is given in Tables Dl to D5. A typical plot of w/W versus x/A 
is shown in Figure Dl. The variance in deflection profile for a con- 
centrated load (C/A = 0.02) and a uniform load (C/A = 1.0) is illus- 
trated in the figure. Other assumed deflection profiles are also 
plotted in order to provide a comparison with other membrane solutions. 

The only reason for assuming a deflected shape of the plate 
is to determine the average strain within the plane of the plate. This 
is accomplished by measuring the length of the deflection profile, 
subtracting the original length of plate, and dividing the difference 
by the original length. If L equals the length of the transverse 
deflection profile, and A is the width of the plate, then the average 
transverse strain is: 

G - L-A (in/in) (D2) 

X A 

In order to determine L from the w/W, x/A relations, the 
deflection profile was divided into a series of straight lines. This 
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A 



Figure Dl. Deflection Profile along the Transverse Axis of a Partially 
Loaded Simply Supported Plate 
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linear approximation of the curve allowed rapid computation of L by 
the computer. In addition, the parameter W/A must be introduced so 
that the strain calculation is related to the actual plate deflection. 
Figure D2 represents a certain segment of the profile curve. . By mul- 
tiplying each segment of A (w/W) by W/A, the arc length A L will 
represent an actual segment of the deflection profile. 




Figure D2 . Element of Deflection Profile' Curve 



Thus 



where 



* L, - ~\/[a (x /A ) ] 2 + [a (w/W) (H/A) ] 2 

= • (5) i + 1 

and 

(U - (5) 1 + 1 

Then, the actual length of the deflection profile will be: 

J 



L = 



A L. 

= i J 



(D3) 



(D4) 



where J equals the number of divisions of the x/A axis. With this 
value of L, the average transverse strain can be found from equation (D2). 

For the case of a square plate, the average strain in the y 
direction will equal the strain along the x axis. If B/A > 1.0, then 

the two strains will be different. The transverse strain will be the 
same as found from equation (D2). The longitudinal strain can be found 
using a similar relation, but the length of the curve must be adjusted 
to reflect the longer dimension of the plate. This can be accomplished 
simply by dividing the W/A term by B/A so that the parameter W/B is 
used in equation (D3) in place of W/A. Then, if is the length of 
the arc along the y-axis, the average longitudinal strain is: 

e y = i/j - B (D5) 

B 



Finally, the principle strain at the mid-point of the plate 
can be determined from equation (CIO) of Appendix C, which is repeated 
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here for convenience. 



6-j = € x + ^ ^ y (D6) 

i - v 2 ' 

This is the value of strain that is presented in Tables Dl to D5 and 
Figure D3 for various values of B/A, C/A, and W/A. The strain is 
referred to as the average transverse strain in the tables because 
only deflections along the x-axis are computed in the analysis. 

Also, it should be noted that Poisson 1 s Ratio is equal to 
0.5 for the calculations. This is the value that is commonly used 
for strains in the plastic range. Both V - 0.3 and 0.5 were used in 
the calculations to test the validity of the membrane formula and 
V = 0.5 generally gave the best results. 
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TRANSVERSE DEFLECTION PROFILE AND. AVG, STRAIN 



B/A= 1,000 POISSON S RAT 10= 0.500 



RATIO 


) OF DE P L . 


COEFFI f 


AT X/A 


TO OFF L 


o COEFF 


IC, 


AT MID-PT 




X/A ,50 


,40 


.30 


. 7 0 




. 10 


C/A 


W (MAX) 














0.020 


0. 0003 5 


1, 0000 


0.91 56 


0, 7268 


0.5034 


0. 


2546 


0. 040 


0. 00141 


]., 9000 


0.9163 


0, 7 2 7 3 


0. 50^0 


0. 


2 549 


0,060 


0. 00316 


V—* 

o 

o 

o 


0.9171 


0. 7293 


0, 50^9 


0. 


? 5 5 6 


0.000 


0, 00559 


I .0000 


0 . 9132 


0 » 7 3 1 2 


0.5061 


0. 


2 564 


0,1 00 


0, 00 36 8 


1 . 0000 


C. 91 96 


9. 7336 


0.5077 


0. 


2 57 4 


0. POO 


0. 03325 


l. 90 00 


0.9281 


0. 7 497 


0.5190 


0. 


7 643 


0.400 


0. 11 500 


loOOOO 


0. 9^31 


0,7828 


0. 559P 


0. 


2328 


0. 400 


0. 20991 


1 , 0000 


0.95Q5 


0. r> 058 


0. 5807 


0. 


3021 


0, 800 


0. ? 833 3 


1. 0900 


0,9537 


0, 3197 


0. 599g . 


0. 


3177 


1 . 000 


0. 3107° 


1 . 0000 


0. 9549 


0.3213 


0.6063 


0 . 


3239 



AVERAGE TRANSVERSE 


STRAIN F' 


OR GIVEN MAX DE 


r L » 




W/A 


0.0P5 


9,950 


0,0 75 


0, 100 


0. 150 


0 


. 200 


C/A 
















0.020 


0.00276 


0,01 102 


' 9. O’ 470 


0 a 0 4 9 6 5 


0, 99668 


0. 


16 3 3 ?. 


o;o4o 


0, 002 76 


0, 01 103 


0.97/7? 


0. 04369 


0,9O 6 7 ? 


0 , 


1 6847 


0.060 


0. 00277 


0. 01 193 


0, 02476 


0.04376 


0. 0 Q 699 


0, 


1 6 369 


0. 080 


0. 00277 


0, oi 1 07 


0 ,02480 


0,04394 


0. 09708 


0. 


16 39? 


0.100 


0, 00278 


0 . C 1 ] 1 0 


0. 02 4.36 


0.04395 


0, 09730 


9, 


16 9 3S 


0. ?00 


0. 90283 


0,01128 


0. 0? 6P3 


0.04467 


0. 09373 


0, 


17186 


0.400 


0. 09295 


0. Cl 175 


0. 02 69? 


0. 04647 


9. 10264 


9. 


17 80.8 


0.60 0 


0. 003 0 5 


0.01217 


9,02723 


0,9 4 804 


0, 19591 


0, 


13 335 


0, 800 


0. 00312 


0. 012-45 


0. 92 79 8 


0. 0491 l 


0. 19810 


0, 


1 3685 


1.000 


0. 0031 5 


0,012 5 5 


0. 0 7 P38 


0,04950 


0. 1089 1 


. 9 


13813 



Table Dl 
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TRANSVERSE DFFl EfTION PR CF 1 1. F AND AVG. STRAIN 



R / A= 1.250 POISSON s RAT I 0 = 0.500 



RATI C 


) OF OFFL 


. C O F f F I C • AT X/A 


TO OF c 1 


I., COEFFIC. 


AT MIp-PT 




X/A =50 


.40 


c 

o 


o 20 




. 10 


C/A 


W (MAX) 














0.070 


0. 00047 


1 . 0000 


0.9715 


0. 7 406 


0.5175 


0 = 


2637 


Oo GAO 


0. 0016° 


1.0000 


0 . 9 7 7 0 


0 . 7 M 4 


0.5130 


0, 


2 640 


0.06 0 


0. 00330 


1. 0000 


0 . 0 7 7 3 


0. 742 7 


0.5138 


0, 


2 64 5 


0. 08 0 


Oo 00673 


lo 00 00 


0. 92? « 


0 » 7 444 


0.519= 


0. 


2652 


0.100 


0. 01 0^.7 


1. 0000 


0.9250 


Oo 7464 


0.5212 


0 = 


2661 


0.700 


0. 04033 


] . ooon 


0. 9??] 


0,7596 


0.5310 


0 = 


2721 


0=400 


0. 14244 


1, 0=00 


0. = 4^6 


Oo 7376 


0 . 5 5 7 8 


0, 


2878 


0.600 


Oo 76325 


1. 0000 


0. 9506 


0. 3065 


0, 5322 


0, 


3 036 


0.800 


0. 37335 


1,0000 


0 o °5?) 


O e - 1 5 6 


0. 5970 


= , 


3156 


1.000 


0. 4M 12 


1 = 000 0 


0,9533 


0. 817Q 


0,6=12 


0. 


3 1=7 



AVFPAGF TRANSVERSE STRAIN for r, I \/ r >\| MAX DEFl . 



W / A 


0.025 


Oo 05 0 


0,0=5 


0, loo 


0.150 


0 


, 200 


C/A 
















0,020 


0, 00246 


0. 0 0=84 


0, 0=205 


0.038=8 


0. 0363= 


0. 


15055 


0,040 


0. 00747 


0. =0=35 


0, 0=206 


0, 03=01 


0, OP 646 


0, 


1 5067 


0. 060 


0. 00=47 


0. 00=86 


0.02210 


0.0 3.30 7 


0, OB 65 8 


0, 


15036 


0.080 


0. =0747 


0, 00=38 


0, =2=14 


0.03=14 


0. 08673 


0. 


15111 


0,100 


0= 00748 


0= 00 99 0 


0,02213 


0. 039=2 


0. 006=1 


0, 


15141 


0.200 


0. 00=5? 


0, 01005 


0. 07 252 


0,0? =81 


0,03816 


0. 


15343 


0.400 


0, 00=61 


0.0104= 


0, 0=334 


0,04124 


0, 0=1 13 


0, 


1 5 84 4. 


0,60 0 


0. =0=69 


0. 0 1 0 7 3 


Oo 0 = 4 = 2 


0 , 0424= 


0. 0=362 


0, 


16240 


0,300 


0. 00774 


Oo 0109 = 


0,92 4 4 ? 


0,04312 


0. 0=512 


0. 


16431 


1.000 


0. 00775 


0.01 ==8 


0, 02456 


0,04335 


0. 0=558 


0, 


16 5 5 5 



Table D2 
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TRANSVERSE DEFLECTION P R TF ILF A N D AVG. STRAIN 
8/ A= 1.500 PH] SSHN S PAT 10 = 0,500 
RATIO OF OFFL, OQ^FFIC. AT X/A TO DEF|_, COEFFIC, AT Ml 





X/A ,so 


,40 


,30 


o 20 


• IC 


C/A 
0, 0?0 


W ( M A X ) 
0. 0°047 


I . 0000 


0,924? 


Oo 7463 


Go 5239 


Oo 2673 


0, 04 0 


0. 001 37 


1. OO 00 


0. °?4 7 


0 » 7 A 7 6 


0. 5?44 


0, 263 1 


0.060 


0. 0 04 1 9 


1.0000 


0. °254 


0, 74 3 3 


0.5251 


0 , 2 6 P 6 


0„ 0 BO 


0. 00741 


1,0000 


0,9263 


0, 75 03 


0, 5261 


0. 26°3 


0. 1 00 


0. 01153 


1 . 0000 


0. °274 


o. 7 5n 


3. 5274 


0. 2701 


0. 200 


0. °4456 


1 o 1000 


0, 9339 


0, 7643 


0, 5364 


0,2756 


0,400 


Oo 1 5 396 


1 « 00 00 


0 . 9 4 5 3 


0,7900 


0, 5609 


0, ?aOO 


0. 600 


0, 30376 


1. 0000 


0, 9 506 


Do 306 3 


0,5829 


0,3042 


O.ROO 


0. 43 710 


l, OOCO 


0, 05 ? a 


0. ° 1 4 7 


0, 5 93 8 


7,3147 


I . 000 


0. 52354 


1,0000 


0, 9 534 


Co 31 66 


Oo 5991 


Oo 3131 



AV F 1 


R AG p TRA’ 


■j S V c R S F 


STRAIN Fnp, GIVE' 


! MAX np-FL, 




W / A 


0.07 5 


0.050 


0.0 75 


0. 100 


0.150 


0 


.700 


C/A 
















0.020 


0. 00730 


0. CC°1 ^ 


0. 02055 


0, O 7 634 


0, OB 054 


0. 


y 403 7 


0.04 0 


0, 00730 


0,00718 


0. 0 2 0 5 7 


0.33637 


0, 03060 


0, 


14Q45 


0,063 


0, 0023 0 


0.0091° 


0 , 0 2 0 a- 0 


0,03642 


C, 03 07 0 


0, 


14062 


0„0 P 0 


0. 00231 


0,00°?] 


0.02063 


0 e 33 64 P 


0. 03 07 4 


0, 


14 0° 4 


0,100 


0, 00231 


0. 7 0 ° 2 3 


Q, 3 2 c 6 8 


0,03456 


0. 08100 


0, 


14111 


0. 200 


0. 00235 


0. 7 3°3 6 


0, 02098 


0, 07703 


0. 07210 


0, 


14294 


0.400 


0. 00243 


7 . 00 96 3 


0, 7? 1 65 


0. 03 331 


0. 03473 


0. 


14725 


0.600 


0. 00249 


0 9 3? QQ4 


0. 02??5 


0. 03 °2 ° 


0. 03678 


0, 


1 5060 


0.30 0 


0. 0° 2 53 


3, 0100° 


0. 0225° 


0.03978 


0, 08300 


0. 


15255 


1 . 000 


0. 00254 


0.01014 


C, 0226° 


0. 04004 


0, 08334 


0. 


15310 



Table D3 



TRANSVERSE D E F L r C T IflNJ P8PFTLF AND AVG. ST 9 A I M 





8/ A = 


1 , 7 c 0 


POI SSOM 


S R AT I 0 = 


= 0,500 




RAT IF 


i OF O c F l 


, CQEFFTC 


AT X/A 


, tq b f f l 


, COFFFIC, AT 




X/A ,50 


. 40 


,30 


. 20 


. 10 


C/A 


VI [ M A v ) 












0.020 


0, 0004 9 


1 , 0000 


0. ° 2 5 5 


o, 7 4 0 R 


0, 5270 


0.2698 


0, 040 


0. 001 06 


1 . 00 00 


0*o260 


0,7506 


0,5275 


0. 2 7 0 1 


0,060 


0, 0044 0 


] , 0000 


0, 9 266 


0. 7517 


0,523? 


0,2706 


0,080 


0*00779 


1.0000 


0*9275 


0, 7 53 2 


0. 5291 


0. 271 2 


0, 1"0 


0. 01212 


1. OPOO 


0. °?85 


0,754° 


0.5303 


C.?720 


0,200 


0. 0460? 


] . 0000 


0 , 0 ?4 8 


0, 7665 


0, 5 3 Q 0 


0. 2772 


0.400 


0, 1 6 R ?0 


1. OQPO 


Go ° 4 5 6 


C. 7910 


0.5624 


0. 2 9 1 0 


0,600 


0. 323' 7 7 


1. 00O0 


0.9507 


0, 8069 


0,5332 


0,3045 


0.800 


0, 4 7 05? 


1 , oooo 


0.9 r ' 2 7 


0. R 1 43 


0.5952 


9. 3143 


1. 000 


0® 671 0? 


1. 00 00 


0,9 53? 


0. 8159 


0.598? 


0,3173 



AVER AG F TRANSVERSE STRAIN F n P GIVEN MAX 0 E F L , 



VI / A 


0. 025 


0.0*0 


0.075 


Oc 1 00 


0,150 


0 


, 200 


C/A 
















0,020 


0,00219 


0. OOP 76 


0. 01963 


0.0347Q 


0, 07638 


0, 


12 394 


0,06 0 


0, 002 9 0 


Go 0 0 8 n 7 


0, 01964 


o. Q7 47 3 


0. 0~ 7 694 


9. 


13404 


0.0 60 


9. 00220 


Go n 0 q 7 8 


0, 01967 


0,03477 


0. 07704 


0, 


1 34 20 


0. 080 


Oo 09? 20 


0o°0 Q 7° 


0. 01970 


0,07463 


0, 07716 


0, 


13 440 


0.100 


0. 00221 


o . c o a a i 


9, 91 974 


0. 03 4 90 


0. 07731 


0, 


] 3465 


0. ?00 


0. 00224 


o , o ^ 


0. 0 2 00? 


0.03538 


0, 07833 


0. 


1 3 6 3 4 


0.400 


0, 00231 


Oo 00 Q ?3 


0, 0 2C6 -7 


0,03652 


0, 03 O' 7 4 


O 


14 030 


0,60 0 


0. 00237 


Oo 00 946 


0, 0? 113 


0.07741 


0, 93240 


0. 


14333 


0,8 00 


0, 00?61 


0. -0^960 


0, O2143 


0,077°? 


0, 98363 


0. 


14506 


1 .000 


0. 0 024 ? 


0,00063 


0. 9? 1 56 


0. 0? 806 


0. 08 396 


0, 


1 4 5 5 1 



Table D4 



D-°T. 



90 



TRANSVERSE D E F L £ C T I n N PROFILE ANC AVG, ST°ATN 





R/A= 


2.000 


P 0 I SSON 


S RATIO- 


‘ 0 , 530 






RATIO OF OE^L. 


C OF n f I r 


. AT X/A 


TO OEF L 


, COEF 


no. 


AT 




X/ A .50 


, 40 


.30 


, 20 




. 10 


C/A 


V ( M A X ) 














0.020 


0. 00050 


1 .0300 


0 . P? 6 ? 


0. 7513 


0, 5285 


0 • 


2 70 3 


0*04 0 


0. 00201 


1 , 0000 


0,9266 


0 . 7 5 3 0 


0. 5290 


0 , 


2 711 


0.060 


0. 00451 


1 a 0000 


0.9273 


0. 7531 


0.5297 


0 , 


2715 


0.080 


0 . 00 BOO 


1 .0300 


0 . 9281 


0 . 75^6 


0, 5306 


0 , 


2 722 


0 a 100 


0 . 012 4 4 


1 . 0090 


0 , 32^1 


0.3563 


0.5318 


0 , 


2729 


0.200 


0. 04 BIB 


1. 0000 


0.9352 


0, 7676 


0, 6 40 2 


r\ 

w O 


2781 


0.400 


0 . 17316 


1 , 9000 


0,9458 


0 . 7916 


0,5631 


0 , 


2 Q 1 6 


0.6 00 


0. 33454 


1 . 0030 


0 . 3507 


0. 8070 


0. 5833 


0 . 


3 047 


0. BOO 


0. 48 864 


1,0000 


0, 9526 


0.814] 


0, 5949 


0 , 


3141 


1 , 000 


0. 596°° 


1 . 0000 


0 . 9 531 


0. 3156 


0,5977 


0 . 


3170 



average transverse STRAIN FOR GIVEN MAX 9 E “ L , 



W / A 


0. 02 5 


0 o 05 G 


0,0 75 


0, TOO 


0. 150 


0 


• ? 0 G 


C/A 
















0.02 0 


0. 00213 


0* 00 548 


0. 01901 


0,03361 


0, 0 746 5 


9. 


1 2*66 


0.040 


0, 00213 


0 a n n 5 4 g 


0. 01 °0? 


0,03363 


3. 07450 


Oo 


12^75 


0,060 


0, 0021. 3 


Oo 0 0 ° 5 0 


0. 01905 


0. 03368 


0, Q7489 


0. 


12*90 


0.06 0 


0.00213 


0 o 00*52 


0. 01 908 


0.03373 


0 , 0-^47 1 


0, 


1 3 0 1 G 


0. 100 


0, 0071 4 


Oo 00° 53 


0,01912 


0.03380 


0. 074B 5 


Oo 


13033 


0.200 


0. 0021 7 


Q a 003 5 5 


C. 01938 


0.034-26 


0, 07533 


0 . 


] *3 194 


0,400 


0, 00224 


0*00893 


0, 02 301 


0,03534 


0, 0781 1 


0 . 


1 3 5 6 a 


0.600 


0 . 002 2 Q 


Oo 0 0015 


0,02049 


0,0361 7 


0. 07986 


0 , 


1 3 8 54 


0.800 


0. 00233 


Go 099?R 


0 . 02077 


0, 03665 


0. 08086 


0 , 


14Q15 


1.000 


0 . 00733 


0 a 0 093! 


0.02084 


0. 03678 


0, 08 ] 12 


0 . 


14056 



Table D5 



Table D6 



LIST OF INPUT VARIABLES FOR COMPUTER PROGRAM FOR. PLATE DEFLECTION 

PROFILE AND TRANSVERSE STRAIN. 



Computer 
Input Variables 


Definition or Symbol in Text 


BAl 


First value of B/A to be used in Eqn. (Dl) 


C2 


Increment of B/A to be added to previous 
value for next set of calculations 


C3 


Last value of B/A desired to be used in 
calculations 


PR 


Poisson's Ratio 


NWA 


Number of W/A values to be used in calculations 


NXA 


Number of x/A values to be displayed in output 


WA 


W/A 


Output Variables 
W (MAX) 


Solution of equation (Dl) for x/A = 0.5, 
the mid-point of the plate 


x/A 


For plate dimensions of Figure 2 


C/A 


For load definition of Figure 1 


W/A 


For deflection definition of Figure 3 



See statements DEFL0004 and DEFL0006 in the program listing 
See Tables Dl - D5. 
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*— « 


CM 


fO 




in 


<5 


X 


00 


O' 


c 




CM 


m 




LO 


o 


X 


CO 


a 


O 
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CM 


CO 




in 


o 


X 


cr 


O' 


o 


H 


o o 


a 


o 


o 


o 


O 


o 


O 


r— H 


rH 


rH 


fH 


•H 


r*-< 






r— i 


#h 


CM 


CM 


CM 


CM 


CM 


CM 


CM 


CM 


CM' 


CM 


CO 


fO 


o o 
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O 


o 


o 
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o o 
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O 


O 


O 
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o 


o 
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F 1= ( { -1 ) **( ( M — 1 ) / 2 ) ) * S I N ( 2 . 0 *GA i v , V A ) /y**5 DE FLOG 3 2 

F2=CGSH( ALPHA-2 *0*GAMMA) ! DEF LOG 33 

F 3= GAMMA# SINH(ALPHA-2.0') : GAMMA) | DEFL003A 

F4= ALPFA*SIMH(2. 0*G A'4M A ) /( 2.0*CSHA) DEFL0035 

Sl=Fl*(1.0-( F2+F3+F4)/CSHA)*SIN (GET A) DEFL0C3 6 
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//11X, »W/A' ,6F8.3,//10X, ’C/AM DEFL0063 

PCA=0 .02 i DEFLC06P 

DC =0 . 02 i DEFLCC70 

DO 75 1=1, 10 DEFL007 1 

WRITER, 103) PCA, (STRN (KK, I ) ,KK= 1,NWA ) DEFL0072 
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A ppendix E 



COMPUTER PROGRAM FOR ICE /PLATE INTERACTION 

The FORTRAN IV computer program to calculate the mid-point 
deflection of a plate and the amount of energy absorbed at the expected 
point of failure is listed at the end of this appendix. The logic of 
the program follows the step-by-step procedure given on pages 37 to 39. 
The values obtained for (C/A) o from step 1 is used as input to the 
program. The applied force F Q and its corresponding W are determined 
in the first part of the program. Then the absorbed energy U is cal- 
culated from the strain relation at the mid-point of the long side of 
the plate. Next, the applied force is increased by a A F and the pro- 
cess is repeated to find the new value of W and U. The output format 
of the program is illustrated in Tables El to E19. The description 
of the input and output variables and their relation to the parameters 
defined on the page of Nomenclature are given in Table E20. 

Similar to a hand calculation of W and U, a constitutive 
relation between the stress and the strain for the plate material is 
required for the computer program. A generalized stress-strain curve 
was used in the program to obtain a stress for a given total strain. 

The general form of the equation is: 

or - d( x + e p ) n * d (e) n (ei) 

where 6p is the plastic strain and 6 is the total strain of the 
material. With the proper choice of coefficients, this equation pro- 
vides a reasonable fit to an actual tensile test stress-strain curve. 



96 



It was found that by setting X equal to the yield strain of a given 

material, and using the boundary condition 0=0 when £ = £ , the 

u u 

exponent n and strength Coefficient D can be determined. For materials 
such as steels which show near constant stress for the smaller values 
of £ , the value used for X should be the strain at the point where 
O' ceases to be constant. For example, the stress-strain curve for 
mild steel exhibits a near constant stress at the yield point up to 
approximately 6 = 0. 02 in/in. Thus X = 0.02 should be used in 
Equation (El). 



Since the generalized expression is always increasing with 
increasing 6 the choice of £ ^ is important to ensure a proper 
curve fit. It is usually difficult to determine from an actual stress- 
strain curve, the actual value of strain that corresponds to @ u * 

The best selection of 6 u is at the point where the O - £ curve first 
starts to level off. The slope of the curve may not actually reach 
zero until just prior to the failure strain. But the ultimate stress 
is usually approached to within a very small error at a strain con- 
siderably smaller than the actual 6 . Therefore, in order to main- 
tain a good fit of the curve between £ = X and £ = 6^, a value 
of £ u should be used that may be less than the actual ultimate strain. 

From the chosen values for X and 6 the values for D 
and n are determined from: 

ln( O Q / O u ) (E2) 

ln( X / 6 u ) 



and 



D = 




(E3) 
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The use of the generalized stress-strain curve allowed the 
incorporation of the average transverse strain program discussed in 
Appendix D into the overall program as a subroutine. The name of this 
subroutine is STRAIN and its listing follows that of the main program 
at the end of this appendix. For an input of the variables B/A, C/A, 
and W/A to STRAIN, a value of 6 j as given by equation (CIO) is cal- 
culated and returned to the main program. Then the^ stress for 6 
is found from the generalized stress-strain function. 



98 



ICE/PLATE INTERACTION 



ICE CRUSHING PRESSURE* 45>? PS I 

PLATE DIMENSIONS (INCHES) 

THICK NESS=0>125 
W I D TH= 54,00 
LENGTH* 84,00 



MATERIAL PROPERTIES 

YIELD STRESS* 56000,0 PS I 
ULT, STRESS* 71200,0 PSI 
YIELD STRA IN=0, 0200 IN/IN 
ULT, ST RAIN* 0,1400 IN/ IN 
POISSON S RATI 0 = 0, 500 

AT f'ID-PT OF LONG SIDE 



F 

(L RS ) 


W 

( IN ) 


S 

( LR/IN ) 


C/A 


STRAIN 

(IN/IN) 


STRESS 

( PSI ) 


U 

( IN-LB) 


205027, 1 


1 , 0068 


7000,0 


o 

o 

o 

r 


0,0230 


56967, 5 


1.309,0 


410027, 1 


3, 8133 


7000, 0. 


u 

o 

o 

o 


0, 0460 


62054, 5 


2851,6 


615027, 1 


5,2830 


7577, 6 


1 ,000 


0 , 066 1 


64903, 2 


4290, 8 


820027, 1 


6, 6620 


801 3.4 


1 , 000 


0,0857 


67012,7 


5741,0 


1025027, 1 


7, 9763 


8366, 2 


1,000 


0, 1 048 


68701, 9 


7201 , 2 



Table El. Plate IB of Reference E3 



ice/PL A 1 E INTERACTION 



ICE CRUSHING PRESSURE- 43.9 PSI 



PLATE DIMENSIONS (INCHES) 
THICKNESS-O. 182 
WIDTH- 54.00 
LENGTH- 84.00 



MATERIAL PROPERTIES 

YIELD STRESS- 41800.0 PSI 
ULT. STRESS 65500.0 PSI 
YIELD STRAIN-0.0175 IN/IN 
ULT. STRAIN-0.1600 IN/IN 
PCISSOM S RAT 10=0. 500 



AT MID-PT OF LONG SIDE 



F 

( LBS) 


W 

(IN) 


S 

( LB/IN) 


C/A 


STRA IN 
( IN/ IN) 


STRESS 
(PSI ) 


U 

( IN-LB ) 


221810.3 


1 . 89 81 


7607. 6 


1 . 000 


0. 0333 


47631.6 


1 586.3 


443810.3 


3.7318 


7742.2 


1.000 


0.0660 


54727.4 


3612. 8 


665810.3 


4. 9 E 42 


8696.5 


1.000 


0.0933 


58702.7 


5474.4 


887810.3 


6.1231 


9439.3 


1 . 900 


0. 1194 


61725.4 


73 71 . 9 


1 1C581 c. : 


7. 1896 


10049 . 3 


1 .000 


0. 1449 


64197.4 


9303 .7 



Table E2. Plate 2A of Reference E3 



ICE/PLATR interaction 



ICE CRUSHING PRESSURE* 241,0 PSI 

PLATE DIMENSIONS (INCHES) 
THTCKNESS=0,244 
VJ I 0 TH = 13.50 

LENGTH* 21,00 



MATERIAL PROPERTIES 

YIELD STRESS* 56000,0 PSI 
ULT, STRE S S= 70700,0 PSI 
YIELD STRA IM=0» 0200 IN/IN 
ULT , STRA IN=0, I 400 IN/ IN 







POISSON S 


RAT I 0 = 0 


, 500 














AT MID 


-PT OF 


LONG SIDE 


F 


W 


S 


C/A 


STRAIN 


STRESS 


U 


( LBS ) 


(IN) 


(LB/ IN) 




( IN/IN) 


(PSI ) 


( IN- LB) 


68323, 5 


0,3255 


I 3664,0 


1.000 


0,0306 


53933,2 


1805, 0 


136523, 5 


0, 6505 


13664,0 


1. 000 


0, 0612 


64028,4 


391 8.7 


204723,5 


0,9714 


13720, 1 


1.000 


0, 0916 


67194, 8 


6153,3 


272923. 5 


1,2257 


1 4496. 1 


1,000 


0.1186 


69307,7 


821«,.7 


341123, 5 


1, 4683 


1512 5, 1 


1,000 


0. I 45 0 


70997, 8 


10294, 8 



Table E3. Plate P2 of Reference E2 
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ICE/PLATE INTERACT I PM 



ICE CRUSHING PRESSURE= 19Q„0 PSI 



PLATE DIMENSIONS (INCHES) 
THICK N ESS = 0 •, 1 19 
W I D TH= 13,50 

LEMGTH= 21,00 



MATERIAL PROPERTIES 

YIELD STRESS^ 56000,0 PST 
ULT, STRESS^ 70700,0 PSI 
YIELD STRA IM = 0» 0200 TN/IN 
ULT, STRA IN=0, 1400 I M/IN 
POISSON S RAT I 0=0, 500 



AT MID-PT OF L n NG SIDE 



F 

(IPS) 


W 

( IN ) 


S 

(LB/IN) 


C/A 


STRAIN 
( I N / 1 N ) 


STRESS 

(PSI I 


U 

( IN-LB) 


53865, 0 


0, 5267 


6664, 0 


1,000 


0,0241 


57278,4 


1383, 1 


107665,0 


1 , 0332 


6783, 9 


1 , 000 


0, 0478 


62161,7 


2971 . 3 


161465,0 


1, 4329 


7336,0 


1 ,000 


0,0688 


64935, 4 


-A469, 3 


215765,0 


In 8 0 ° 1 


7746, 4 


1 > 000 


0,0897 


66990, 1 


5980,3 


769065 , 0 


2, 1 707 


8069, 4 


1 , 000 


0, 1094 


63633, 4 


7506,0 



Table E4. Plate P3 of Reference E2 



Tf E/PLATE INTERACTION 



ICE CRUSHING PRFSSURE= 120,, 0 PSI 



PLATE DIMENSIONS (INCHES) 
THICKNFSS=0,105 
WIDTH= 13,50 
LENGTH 3 21,00 



MATERIAL PROPERTIES 

YIELD STRE SS= 32000,0 PSI 
ULT, STRESS= 46100,0 PSI 
YIELD STRA IN=0. 0200 IN/IN 
ULT, STRA TN=0„ 1 600 IN/IN 
POISSON S RATIO -0,500 

AT MID-PT OE LONG SIDE 



F 

( LES) 


W 

(IN) 


S 

(LB/ IN) 


C/A 


STRAIN 
( I N / I N ) 


STRESS 
(PSI ) 


U 

( IN-LB) 


34020,0 


0, 6592 


3360,0 


1 ,000 


0> 0267 


33662 , 7 


898,4 


68020,0 


1, 2120 


3 65 3, 6 


1 , 000 


0,0511 


37724, 1 


1926 o 3 


102020,0 


1 , 6380 


4054, B 


1,000 


0, 0727 


40136, 2 


2917, 2 


136020,0 


2,0307 


4360, 6 


1 , 000 


0 0 ° 3 6 


419 57, 1 


3 9 2 6 ? 3 


170020,0 


2,4000 


4611,8 


1 >000 


0,1140 


43434, 5 


4950,2 



Table E5. Plate P4 of Reference E2 



ice/plate interaction 



ICC CRUSHING PR E SS UR E = 87,0 PSI 



PLATE DIMENSIONS ( INCHES ) 
THICK NESS=Cj068 
W I OTH= 13, >50 

LENGTH* 21,00 



MATERIAL PROPERTIES 

YIELD STRESS* 32000,0 PSI 
ULT, STRESS* 41800,0 PSI 
YIELD STRA IN = 0. 0200 IN/IN 
ULT, STRA IN=0, 1600 IN/IN 







POISSON S 


RATI 0=0 


,500 














AT MID 


i-n>T OF i 


LONG SIDE 


F 


W 


S 


C/A 


STRAIN 


STRESS 


U 


(IRS) 


(IN) 


(LR/ IN) 


• 


( I N / I N ) 


(PSI ) 


( IN-LR) 


24664,5 


0, 73 79 


2176,0 


1,000 


0,0193 


31364, 3 


616, 6 


492 84 > 5 


1,3536 


2370, 4 


1 , 000 


0,0370 


34627,5 


1280,0 


73904,5 


1, 8709 


2571 , 6 


1,000 


0,0532 


36285, 5 


1930, 3 


98524, 5 


2,3545 


2724, 2 


1 , 000 


0,0690 


37515, 4 


2536,8 


12-3144, 5 


2, 8160 


2846, 9 


1 ,000 


0, 0P44 


33502,4 


3240,6 



Table E6. Plate P6 of Reference E2 



ICE/PLATE INTERACTION 



ICE CRUSHING PRESSURE* 184,2 PSI 

PLATE DIMENSIONS (INCHES) 
THICKNESS=0,119 
WIDTH= 13,50 
LENGTH* 21,00 



MATERIAL PROPERTIES 

YIELD ST R F SS = 42000, 0 PSI 
ULT, STRESS- 61000,0 PSI 
YIELD STRAIN*©. 0200 IN/IN 
ULT , STRAIN =0,1 6 00 IN/IN 
POISSON S RATI 0=0, 500 

AT MID-PT QP LONG SIDE 



E 

UBS) 


W 

( I M ) 


S 

( LB/ IN) 


C/A 


STRAIN 

(IN/IN) 


STRESS 

(PSI) 


U 

( IN-LB) 


52220)7 


0, 6802 


4998, 0 


1,000 


0,0312 


454°2 , 8 


1420,0 


104420) 7 


l, 2409 


5478, 4 


1 , 000 


0,0505 


51074, 1 


3037,9 


156620,7 


1 , 6747 


6088,3 


1 ,000 


0, 0846 


5441 0, 7 


4604, 6 


208820, 7 


2, 072 6 


6559, 1 


1,000 


0,1 030 


56928, 9 


6198, Q 


261020,7 


2, 4442 


6952, 4 


1.0 00 


0, 1325 


58969, 9 


7813,6 



Table E7. Plate P8 of Reference E2 



ice/plate interaction 



ICE CRUSHING PRESSURE = 100.0 n S I 



PL AT F DIMENSIONS (INCHES) 
TH I C KMES S = 0 f, 188 
WIDTH= 13,00 
UENGTH= 36,00 



MATERIAL PROPERTIES 

YIELD STRESS^ 38500,0 PSI 
ULTo STRESS^ 59000,0 PSI 
YIELD ST RA I N=0, 02 00 IN/ IN 
Ul.To STRA IM = 0, 1600 IN/TN 
POISSON S RAT 10 = 0,500 



AT MID-PT OF LONG SIDE 



F 

(IPS) 


W 

( IN ) 


S 

< LG/ INI 


C/A 


STRAIN 

(IN/TN) 


STRESS 
( PS I ) 


U 

( IN-LB) 


46800, 0 


0,2843 


7238,0 


1,000 


0, 0205 


38691 , 0 


792,7 


1 40300,0 


0, 8524 


7238, 0 


9 

o 

o 

o 


0,0614 


48472,2 


2977, 2 


233800,0 


1, 2766 


80 53, 7 


1 , 000 


0, 0968 


53216, 3 


5151,0 


327300,0 


1 , 622^ 


8869,0 


1,000 


0, 1292 


56462, 5 


7292, 5 


42.0800, 0 


I , 9440 


9519,0 


1 » 000 


0, 1605 


59041 , 5 


947 9 , 1 



Table E8. Plate no 



2 of Reference El 



ICE/PLATE? INTERACTION 



IGF CRUSHING D RFSSU^r= 100,0 PS I 

PLATE DIMENSIONS (INCHFS) 

TH 1 CKNFS S = 0-, 18 p 
WI D TH= 18,00 

LFMGTH= 19,28 



MATERIAL PROPERTIES 

YIELD STRF SS= 28800,0 PST 
ULT. STRFSS= 89000,0 PSI 
YIELD ST P A I N=0 , 0? 00 IN/IN 
ULT,. STRAIN =0. 1600 IN/ IN 
POISSON S RAT I 0 = 0, 800 

AT MID-PT OF LONG MQE 



F 

(LBS) 


W 

( IN ) 


S 

(LB / IN 1 


C/A 


STRAIN 
( I N / 1 N ) 


stress 

( PSI ) 


IJ 

( IN-LB) 


25025,0 


0, ??65 


72 3 8a 0 


1 a 0 0 0 


0.0178 


37590, A 


660,2 


58825, 0 


Oo 555 8 


72 38., 0 


1 a 000 


0 , 04 1 8 


44799, 8 


18/4,6 


92625,0 


0,8752 


7238,0 


1,000 


0 , 065 c 


49175,3 


3260, 1 


1 264? 5, C 


1, 1059 


781 8, 5 


1 a 0 0 0 


0, 086 4 


51983,4 


44Mf>~8 


160225,0 


1 , 308° 


8371,6 


la 000 


0,1057 


54186 , 4 


57 2 ! a 4 


1 9402 5, 0 


1, 5001 


8 345a 6 


1 a 000 


0 , I 246 


56044- 7 


6981 ,3 


227825,0 


1 , 6 8 ? 9 


9262,3 


1,000 


0,1431 


57660, 8 


8260, 5 


261695,0 


1 , 8569 


96 35, 4 


1 , oon 


0,-1613 


59096 , 2 


95 U) „ 8 


295425,0 


2 , 02 5 6 


9974, 1 


la 000 


0. 1792 


60391, 2 


1 087 6 , 3 



Table E9. Plate no. 3 of Refei.-nce El 



ICE /PL A IF ! NTf RAC HUM 



ICE CRUSHING PRESSURE- 100.3 
PLATE C I I EM S I G NS (INCHES) 



1 H ICKNE SS=0 . I BS 
■'! I OTH = 13.00 

LENGTH^ 13/)0 



PS I 



MATERIAL PROPERTIES 

YIELD STRESS^ 33 3 00 . C PS I 
UL r . STR ESS- 59000. 0 PSI 
YIELD ST RAIN= 0.0290 1M/IM 
ULT. ST PA I i\= 0.1 6 00 IN/ IN 
PR IS SON S RAT in = 0.500 

AT M IO-PT PF LONG SIDE 



F 


W 


s 


C / A 


STRA I N 


STRESS 


U 


(l BS ) 


( IM ) 


( LB/ IN ) 




{ I N / I N ) 


(PSI) 


( IN-LB) 


165 CC. 0 


0. 17 17 


72 3 3. 0 


1 . COO 


0. 0142 


3 53 7 C . 3 


508 .3 


50703.0 


0.5151 


72 3 3 . C 


1 . 003 


0. C425 


44544.9 


1 9 1 0 . 5 


845CC. C 


0. 3 A 8 2 


732 5 . 9 


1 . 000 


0.0704 


4 9 8 4 R . 0 


3 5 C 3 . 9 


1 13300.0 


1 . C 7 14 


3 0-3 1 . 6 


1 . 000 


0. 0936 


52355.1 


4 9 4 9_ . 2 


1-52 1 C C . 0 


1 .23 3 5 


36 30. 3 


1 . 000 


0.1 162 


55247.1 


64 17. 9 


1 35003,0 


1 . 'if (6 


9155.1 


l . 900 


0 .133 1 


57 24 1 .9 


7903 .8 


2 19 7 30 . C 


1.57 50 


5665.2 


1 . 3 30 


0. 1 5 1 '6 


58966. ? 


9 4 C 3 . 3 


2 53 50 C. C 


1 . 0555 


1 0346.3 


1.000 


0. 13 07 


6 049 1 . 2 


1C530.C 


2 8 7 3 C C . 0 


2. 0296 


10405 . G 


1 . ccc 


C. 2 Cl 6 


6 1362.6 


12467.5 



Table E10. Plate no. 5 of Reference El 



ICE/ PL ATE INTERACTION 

ICE CRUSHING PRESSURE^ 100,0 PSI 

PLATE DIMENSIONS (TNCHFS) 
THICKNESS=0,250 
WIOTH= 13,00 
LENGTH^ 26,00 



MATERIAL PROPERTIES 

YIFLD STRESS- 66000,0 PSI 
ULT. STRFS S= 83000,0 PSI 
YIELD STRA !N=0, 0200 IN/IN 
ULT, STRAIN=0ol200 IN/ IN 
POISSON S RAT I 0 = 0 , 500 

AT MID-PT OF LONG SIDE 



F 

( L PS) 


W 

( IN ) 


s 

( LB/TN) 


C/A 


STRAIN 
( I N/IN) 


STRESS 

(PSI) 


U 

(IN-LB) 


33800,0 


0,1167 


16500, 0 


1,000 


0, 0114 


61412, 9 


699,4 


1 01 300,0 


0, 3497 


16500, 0 


1,000 


0,0341 


70669,7 


2412,0 


16HP00 , 0 


0,582? 


16500, 0 


1 , 000 


0, 05 6° 


75439, 6 


4290, 5 


236300, 0 


0., 81 57 


16500,0 


1,000 


0, 0706 


78756,4 


6270,2 


303 800, 0 


1 ,0289 


16817, 1 


1,000 


0,1013 


81223, 2 


8223,8 


371300,0 


1,2075 


17513, 8 


1 ,000 


0,1212 


83107, 2 


10074, 1 



Table Ell. Plate no. 4 of Reference El 
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TCE/PIATF INTERACTION 



F 

(I PS) 

42 3 3 6., 0 

1 42336.9 

242335. 9 

342335. 9 

442335.9 

542335.9 

642335.9 



ICE CRUSHING PRE SSUR E= 600,0 PSI 



PI ATE DIMENSIONS (INCHES) 
THICK NESS=0 J 750 
W I DTH= 24,00 
LENGTH* 48,00 



MATERIAL PROPERTIES 

YIELD STRESS* 38590,0 PSI 
Ul.T, STRFSS* 59000, 0 PSI 
YIELD ST PA IN=0 ,0200 IN/IN 
ULT, STRA IN=0, 1600 IN/IN 





POISSON S 


RATI 0=( 


)„ 500 












AT MID 


i-PT OF 1 


LONG SIDE 


w 


s 


C/A 


STRAIN 


STRESS 


U 


(IN I 


(LB/ IN ) 




( IN/IN) 


(PSI) 


( IN-LB) 


0, 8001 


2 3 87 5, 0 


0,350 


0 , O 7 5 7 


5 0 5 9 5 > 7 


3829,4 


1* 2013 


23875, 0 


0, 642 


0, 1388 


57301 , 5 


7952,3 


1,3350 


28575, 0 


0,83 7 


0,1511 


60518,4 


10953, 8 


1, 36 36 


28875, 0 


0, 995 


0,215? 


62702 , 8 


13495,2 


1, 5424 


28875, 0 


1,000 


0, 2446 


64374, 1 


15749, 1 


U 7079 


28875, 0 


1,000 


0,2709 


6573 5 . 0 


17807,3 


1, 830 8 


29716, 0 


1,000 


0, 2Q?a 


66775, 7 


19527,6 



Table E12. Arbitrary Plate no. 1 
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TCE/PL4TE INTERACTION 



ICE CRUSHING PPFSSURE= 600,0 PS I 

PLATE DIMENSIONS (INCHES) 
THICKNESS-1,000 
WIDTH* 24,00 
L ENGT H= 48,00 



MATERIAL PROPERTIES 

YIELD ST R F SS= 38500,0 n SI 
ULT„ STRESS^ 59000,0 PS I 
YIELD STRA IN = 0,0200 I M/IN 
IJLT, STRAIN=0.1 600 IM/IN 







POISSON S 


RAT I 0=0 


, 500 














AT MID 


i-PT OF 


LONG SIDE 


F 


W 


S 


C/A 


STRAIN 


STRESS 


U 


( LBS) 


(IN) 


( L R / I N ) 




( I M / IN) 


( PS I) 


( IN-LB) 


104543,9 


0, 8259 


33500,0 


0,550 


0,1189 


55515,0 


6602, 8 


204543, ° 


0, 9 7 5 6 


38500, 0 


0, 769 


0, 1664 


5 q 4' 7 4 , 3 


9894,4 


304543 , Q 


1.021 J 


38500, 0 


0,939 


0. 2030 


61954 , 5 


12576,6 


40454? , 9 


1„ 1063 


35500, 0 


1 „ 000 


0,2340 


63786, ^ 


149 2 3,-8 


504543, 9 


1.2355 


38500, 0 


1 , 000 


0, 261? 


65249, 5 


1 7 048 o 8 


604543,9 


1,3524 


33500,0 


I. ,000 


0, 2860 


66471 , 7 


19011, 6 



Table E13. Arbitrary Plate no. 2 



JCF/PLATF: tntfpaction 



ice crushing pressure* 6oo«o rsi 



PLATE DIMENSIONS (INCHES) 
THICKNESS=1 .?50 
W I DTH= 2 4,00 

LENGTH* 4 8, 00 



MATERIAL PROPERTIES 

YIELD ST RE SS= 38500,0 PSI 
ULT. STRESS* 59000,0 PSI 
YIELD STRA IN =0,0200 IN/ IN 
ULT, STRA IM = 0. 1600 IN/ IN 
POISSON S RAT I 0=0 , 500 



AT MID-PT OF LONG SIDE 



F 

( L R S ) 


W 

( IN ) 


S 

( 1. B / I N ) 


C/A 


STRAIN 
( IN/IN) 


STRESS 
( P S l ) 


U 

( IN-LR) 


279935 , 7 


0, 812° 


48125,0 


0, POO 


0, 1946 


61421,0 


11 954., 0 


379935, 7 


0,8577 


43125, 0 


1 , 000 


0, 2267 


63377, 1 


14369., 9 


479935,7 


0, 9640 


48125,0 


I ,000 


0, 2548 


64 9 1 5 , 5 


16542,7 


5 7 9935, 7 


lo 05 ot 


4°1 2 5, 0 


I , 000 


Oo 2301 


66188,8 


1 8541,4 


679935,7 


1, 1 474 


48125,0 


1 , 000 


0, 3033 


67 27 8.. 5 


20406 , 9 



Table E14. Arbitrary Plate no. 3 
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ICE/PLATE INTERACTION 



ICE CRUSHING PRF SSURF.= 600,0 PS I 



PLATE DIMENSIONS (INCHES) 
THICK NFSS=1 ,500 
W I DTH= 24» 00 
L ENGTH= 4 6,00 



MATERIAL PROPERTIES 

YIELD STRESS^ 38500,0 PSI 
ULT, STRESS^ 59000,0 PSI 
YIELD STRA IN = 0,0200 IN/IN 
ULT, STR4IN=0. 1600 IN/IN 
POISSON S R A T I Q = 0 , 500 



AT MID-PT OF LONG SIDE 



E 

(LBS) 


W 

( IN ) 


S 

( 1 B/ IN ) 


C / A 


STRAIN 

(IN/IN) 


STRESS 

(PSI) 


IJ 

( IN- LB) 


691200,0 


0, 681 7 


5 7 750 , 0 


1,000 


0,2162 


62 764, 0 


185 7 ?,6 


701200,0 


1,0314 


57750, 0 


1,000 


0, 3272 


63333, 2 


22358,5 


891200,0 


1 -> 09 4 7 


57750, 0 


1,000 


0,3473 


69173, 1 


24021, 1 


991200,0 


1, 1545 


57750, 0 


1 > 000 


0, 3662 


69932,4 


25611,0 



Table E15. Arbitrary Plate no. 4 



TC F/ PI. AT E INTERACT I ON 



ICE CRUSHING PRRSStJRF* 600,0 PS I 

PLATE DIMENSIONS (INCHES) 

THICK N ES S= 0,750 
WIOTH= 24,00 
LENGTH^ 48,00 



MATERIAL PROPERTIES 

YIELD STR ESS= 66000,0 PSI 
ULT. STRF S S~ 83 000,0 PSI 
Y I C LD STRA IN= 0, 02 00 IN/ IN 
Ul_T, STRAIN =0,1 200 IN/IN 
POTSSON S RAT I 0=0, 500 

AT MID-PT OF LONG SIDE 



F 

(LBS ) 


W 

( IN) 


S 

( LB / IN ) 


C/A 


STRAIN 
{ IN/ IN) 


STRESS 

( PSI ) 


U 

(IN-LB) 


104543, 9 


0, 6423 


49500, 0 


0, 550 


0, 0694 


77382 , 3 


5368, 0 


254543,9 


0,^833 


49500, 0 


0,858 


0,1003 


81914,0 


8867,9 


404543,9 


0,8605 


49500, 0 


1,000 


0, 1365 


84377,4 


11515, 8 


554543, 9 


1, 0074 


49500, 0 


1,000 


0,1598 


8 60 c 6 , 6 


13757,4 


704543,9 


1,1355 


49500, 0 


1 , 0 00 


0,1801 


07425, 0 


15746, 1 


854543, 9 


1, 2506 


49500, 0 


1 ,000 


0,1984 


88510,9 


17556,9 



Table E16. Arbitrary Plate no. 5 
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ICE/PLATE INTERACTION 



F 

( LCS ) 

221183.9 

371183.9 
5211 £3. 5 

671183.9 
.821 1 83.9 

971183.9 



ICE CRUSHING PRESSURE^ 600 .0 PS I 



PLATE DIMENSIONS (INCHES) 
THICK NESS = 1 .000 
W ID TH= 2 A « 00 
LENGT H= A- 8 . 0 0 



MATERIAL PROPERTIES 

YIELD STRE S S= 66000.0 PSI 
ULT. STRE S S= 83000.0 PSI 
YIELD STRA IN=0.C200 IN/IN 
ULT. STRA IN=0. 1200 IN/IN 





POISSON S 


RAT 10 = 0 


w 


S 


C /A 


( IN ) 


( LB/ IN) 




0. 5765 


66 COO. 0 


0. 800 


0.6182 


66000.0 


1.000 


0.7325 


66000. 0 


1.000 


0.8313 


66000. 0 


1. 000 


0.9195 


66000 .0 


1.000 


C. 9999 


66000. 0 


1,000 



.500 

AT MIO-PT OF LONG S IDE 



STRAIN 
( I N / I N ) 


STRESS 

(PSI) 


U 

(IN-LB) 


0.1009 


81181.3 


8192.5 


0. 1307 


83914.2 


1 C 970. 2 


0.1549 


85755.7 


13284.4 


0. 1758 


87154.? 


15321 .2 


0. 1944 


Be 28 5. 8 


17167. C 


0.2115 


89238.2 


18370.6 



Table E17. Arbitrary Plate no. 6 
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ICE/PLATE INTERACTION 



ICE CRUSHING PRE SSURE= 600,0 PS I 

PLATE OT MENS IONS (INCHES) 

THICK NESS=1 >250 
WIDTH* 24,00 
LENGTH* 40,00 



MATERIAL PROPERTIES 

YIELD STRF SS= 66000,0 PS I 
ULT, STR ES S= 8 3000, 0 PS I 
YIELD STRA IN=0,0200 IN/IN 
ULT, STRA IM = 0, 1 200 I N/ I N 
POISSON S RAT I 0=0, 500 

AT MID-PT OF LONG SIDE 



F 

(LPS) 


W 

(IN) 


s 

(LR/ IN ) 


C/A 


STR A I N 
( IN/I Ml 


STRESS 

(PSI ) 


U 

( IN-LB ) 


601 2 0 0 r, 0 


0,4772 


82500,0 


1 ,000 


0,1261 


83531, 9 


10537, 1 


891200,0 


0, 7663 


82 500,0 


1,000 


0, 2026 


88740 , 0 


17977,7 


1 001 2 00,0 


0 , 8470 


82500,0 


1, 000 


0,2241 


89905, 7 


20152,2 


1201200,0 


0 , °? 2 4 


82500, 0 


1 , 000 


0, 2438 


00878, 6 


22158, 6 


1491 200.0 


1, 020 5 


82500, 0 


1 ,000 


0.2721 


92164,9 


25082,4 



Table E18. Arbitrary Plate no. 7 
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ICE/PLATE INTERACTION 



ICE CRUSHING PRESSURE^ 600.0 PSI 

PLATE DIMENSIONS (INCHES) 

THI CKNESS=1 , 500 
W I DTH= 26,00 



LENGTH= 40,00 



MATERIAL PROPERTIES 

YIELD S T RE SS= 66000,0 PSI 
ULT„ S TR E S S= 33000,0 RSI 
YIELD STRA IN=0, 0200 IN/IN 
ULT, STRA IN=0, 1200 IN/IN 
POISSON S RAT 1 9 = 0, 500 

AT MID-PT OF LONG SIDE 



F 

( L P S ) 


W 

( IN ) 


S 

(14/ IN ) 


C/A 


STRAIN 

( I M/IN) 


STRESS 

(PSI) 


U 

( IN-LB) 


601200 , 0 


0,3977 


99000,0 


1,000 


0,1 261 


83531 , 9 


1 0 5 3” 7 , 1 


8 91 2 00,0 


0,6386 


99000,0 


1,000 


0, 2026 


38740, o 
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Table E19. Arbitrary Plate no. 8 



Table E20 



LIST OF INPUT VARIABLES FOR THE COMPUTER PROGRAM 
FOR PLATE DEFLECTION /LOAD RELATION 
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(See Tables El to E19) 



All of the variables displayed in the output are self explan- 
atory or of the same notation as used in the text. The STRAIN at the 
mid-point of the long side is equal to £ t in equation (19) of the 
text. The STRESS at the same point is equal to 



See statements FAIL0003 and FAIL0005 in the program listing. 
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GC TC 5 7 STRN003? 

STN= ( ARC- . 50 )/ .50 ! STRN00R3 

YSTN={YARC-.5C)/ .50 1 STRN0034 

STRN=(STN + PR*YSTN) / ( 1 . 0-PRY*2 ) • STRN0035 
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